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Abstract 

We use a mix of analytic and numerical methods to exhaustively study a class of 
asymptotically global AdS solitons and hairy black hole solutions in negative cosmological 
constant Einstein Maxwell gravity coupled to a charged massless scalar field. Our results 
depend sensitively on the charge e of the scalar field. The solitonic branch of solutions we 
study hits the Chandrashekhar limit at finite mass at small e, but extends to arbitrarily 
large mass at larger e. At low values of e no hairy black holes exist. At intermediate 
values of e hairy black holes exist above a critical charge. At large e hairy black holes 
exist at all values of the charge. The lowest mass hairy black hole is a smooth zero 
entropy soliton at small charge, but a (probably) singular nonzero entropy hairy black 
hole at larger charge. In a phase diagram of solutions, the hairy black holes merge with 
the familiar Reissner-Nordstrom— AdS black holes along a curve that is determined by the 
onset of the superradiant instability in the latter family. 
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1 Introduction 

The AdS/CFT correspondence maps asymptotically AdS solutions of Einstein's equations to 
states of a dual conformal field theory ^-[3|. Static bulk solutions map to 'phases' of the CFT. 
A thorough investigation of all static gravitational solutions consequently permits complete 
understanding of the phase structure of the dual large N field theory, and so is of considerable 
interest. 

In this paper we follow |4] to study AdS Einstein Maxwell gravity interacting with a charged 
massless scalar field. The bulk theory we study is governed by the Lagrangian 
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where = — ieAfj_, e, the charge of the scalar field, is a free parameter, and we have 



set the cosmological length i = 1. The action (1.1), sometimes called the AdS Abelian Higgs 



model, may be regarded as a simple toy model for the charged scalar dynamics of systems 
that appear in concrete examples of the AdS/CFT correspondence (see e.g. [5]-[7]). There 
has been great interest in this model since it was realized that it allows for a phase transition 
between the familiar Reissner-Nordstrom— AdS (RN AdS) black holes and hairy black holes, i.e. 
solutions with a scalar condensate [8]. In the AdS/CFT dual to this system, such condensate 
spontaneously breaks the global U{1) boundary symmetry, and so corresponds to a superfiuid 
phase of the boundary CFT. Bulk solutions that asymptote to AdS with planar sections (i.e. 
Poincare patch AdS) are dual to phases of the dual field theory on spatial R^, and have been 
extensively explored within the context of the AdS/CFT correspondence to learn something 
about condensed matter phenomena [9, 10 ; see [ll}jl3] for reviews. The solutions we study 



in this paper are dual to superfiuid phases of the boundary field theory on a spatial 5"^. The 
results of this paper must, of course, reduce to the planar model in the limit that the radius of 
the S^ is taken to infinity; we will see below that this is indeed the case. 

In this paper we perform a thorough analysis of the properties of three classes of static 
charged solutions of the action ( |1.1 ). The solutions we study all asymptote to global AdS^ 



space (i.e. AdS^ with spherical sections), and so map to phases of the 'dual large N theory' on 
S^. We investigate the properties of our solutions as a function of their mass, charge, and also 
of the free parameter e in our model. 

The first class of solutions we study in this paper is the set of 'ground state solitons' 
discovered in [ij. These 'solitons' are static lumps of a scalar condensate in global AdS^ space. 

At any given value of e these solitonic solutions appear in a one parameter family labeled 
by their charge; the mass of solitonic solutions is determined as a function of their charge. At 
infinitesimal values of the charge these solitonic solutions are extremely simple; they reduce 
to AdS^ space perturbed by the lowest energy linearized mode of the scalar fieldj^ In [i] this 
linearized solitonic solution was used as the starting point of a perturbative construction for 
solitonic solutions in a power series in their charge Q. In this paper we have used numerical 
techniques to extend the construction of [2] to arbitrarily large values of the charge Q. Using 
the analytic construction of [4] as the starting point of our numerical construction, we slowly 
iterate to larger values of the charge. Our results agree perfectly with those presented in |4] at 
small charge. However, we find some surprises at large charge. 

Our numerics, presented in section [3| show that the quahtative properties of solitonic solu- 
tions differ depending on the value of the parameter e. For e > Csoicrit, where Csoicrit is a critical 
charge, the solitonic branch of solutions continues all the way to infinite charge; in other words 
there exist solitons at arbitrarily large charge. On the other hand for e < esoicrit, the solitonic 
branch of solutions terminates in a 'Chandrashekhar' singularity at a finite value of the mass 
Merit and charge Qcrit (see [T] for similar behaviour in a closely related context). The approach 
of the solitonic solution to the critical value is characterized by a spiral behaviour and critical 
exponents familiar from the study of boson stars (see (?]). Numerically we find that egoicru is 

very close to a value that we will encounter again below. 

As e is cranked up to esoicrit from below, our numerical results appear to indicate that M^rit 
and Qcrit diverge in a manner proportional to ln(e — Csoicrit)- We have drawn this conclusion 
by noting that a graph of Qcrit versus Ini^esoicrit — e) is very nearly a straight line at small 
Gsoicrit — e. Because of the slow growth of the logarithm, however, we have not been able to 



^The solitonic solutions map to 'bose condensates' in the dual CFT on S^. 

^The harmonic time dependence of the mode is gauged away yielding a static solution, at the expense of 
turning on the value of Aq at the boundary. 
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verify this divergence in a more straightforward manner by simply checking that Qcrit grows 
arbitrarily large as e is taken arbitrarily close to egoicrit- In the absence of analytic understanding 
here, our results on this point should be treated as tentative, subject to further confirmation; 
in particular we have not really ruled out the (unlikely sounding) possibility that Qcrit remains 
finite at e = esoicrit- 

The second type of solutions we consider in this paper is the very well known RN AdS 
charged black hole solutions. We study the stability of these solutions to the condensation of 



the scalar field in (1.1), as a function of the black hole mass, black hole charge, and parameter 
e. Our analysis of this question (which uses a mix of analytic and numerical methods) yields 
the following results. RN AdS black holes are never unstable to scalar condensation when 
< 3. On the other hand when 3 < < y, near-extremal black holes display an instability 
at large enough charge, i.e. when Q > Qo{e'^). In the range 3 < < y the minimum 
charge for instability, (5o(e^), is a monotonically decreasing function of e^, with Qo{3) = oo and 
Qo(y) = 0. Black holes with Q > Qo(e^) - and for all values of the charge when > y - are 
unstable in a mass band around extremahty. The function Qo{e^) itself appears to undergo a 
phase transition at = y in this range, associated with a change in the qualitative nature of 
the instability mode at this value of e^. We pause momentarily to discuss this in more detail. 

The instabilities that leads to scalar condensation about RN AdS black holes are easiest to 
study in the case of extremal black holes and occur in at least two varieties. The first kind of 
instability is localized entirely within the near horizon 74^5*2 geometry of the extremal black 
holes. This instability occurs because the scalar field acquires a near horizon effective mass that 



violates the Breitenlohner-Freedman bound of the near horizon AdS2 14 , and turns out to be 
the dominant instability for extremal black holes whose horizon radius is larger than the radius 
of curvature of the ambient AdS^ [15j. However, for solutions asymptoting to global AdS (as 
opposed to planar solutions) the near horizon region is not always the determining factor for 
instability. Ref. [Jj] found that for small black holes (compared to the AdS scale) the instability 
mechanism is not dominantly near horizon. The instability is better thought of as an AdS 
superradiant instability of charged scalar fields. ^j^ As we argue below, for extremal black holes, 
the switch between these two different types of instabilities leads to a non analyticity in the 
function Qo{e^) at a particular finite value of e^, most likely = y. We mention that there 
is a closely related ongoing programme to find the phase diagram associated to the AdS scalar 



superradiance for rotating black holes, rather than charged 20 26 . In the future it would be 
interesting to combine these studies and eventually find charged rotating AdS hairy black holes. 

Consider an unstable RN AdS black hole of charge Q and mass M. Triggering the insta- 
bility leads to a decay process that presumably settles down, at infinite time, to a stable static 
configuration of mass M and charge Q. By the Hawking area increase theorem the resultant 
configuration has an event horizon and so is a black hole. As the instability involves conden- 
sation of the scalar field, the resultant black hole solution is immersed in a sea of scalar hair. 
Such solutions are sometimes called 'hairy' black holes, and must exist at every value of M 
and Q for which RN AdS black holes are unstable. Small hairy black holes in the Lagrangian 



(1.1) were constructed perturbatively (in an expansion in the charge of the solutions) in |4]. 
The key observation in [2] is that infinitesimally small hairy black holes may accurately be 
thought of as a 'superposition' of a small RN AdS black hole and the scalar soliton. Using this 
superposition configuration as a starting point, hairy black hole solutions may be constructed 



Recall that a mode e of a scalar field with charge e increases its amplitude by scattering off a charged 



black hole with chemical potential ^ if w < e/Lt 16 19. In an asymptotically global AdS spacetime, this leads to 
an instability since the outgoing wave is reflected back onto the black hole and scatters again, further increasing 
its amplitude [4j[20j[21 . 
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in a perturbative expansion in the charge of the black hole solutions. Working in the small 
charge limit at any fixed value of > y, the authors of j4j found that the lowest mass black 
hole solutions are simply the smooth, zero entropy solitonic solutions described earlier in this 
introduction. These configurations may be thought of as the infinite temperature limit of hairy 
black holes. 

One of the main results of the current paper is that the nature of hairy black holes, in 
the zero temperature limit, changes discontinuously as a function of the charge (see Section |6 



when > y. At low values of the charge - in fact for all Q < Qc2{.^^)^ where Qc2(e^) is a 
monotonically increasing function with (5c2(^) = - the lowest mass hairy black holes indeed 
reduce to the solitons, in agreement with the perturbative results of j4]. For Q > (5c2(e^), 
however, the lowest mass hairy black holes do not reduce to the soliton. The new extremal 
solutions at these charges appear to have finite entropy and zero temperature. They may well 
be singular in the strict extremal limit; in fact, regular extremal hairy black holes are not 



allowed in the system under study 27 . It turns out that for = y (1 + 6*), Qc2(e^) ~ C'(^). 
As a consequence, the phenomenon described in this paragraph - the discontinuous change in 
the nature of extremal black hole solutions as a function of their charge - is visible at small 
charges in a perturbative expansion in 9 (see Section |6]). At larger values of scalar charge e, 
this phenomenon is also visible in our numerical construction of these hairy black holes (see 
Section l6l). 



We now turn to the spectrum of hairy black hole solutions for 3 < < y. As we have 
explained above, in this range hairy black hole solutions exist only for Q > Qole"^). When 

= y(l — 9), it turns out that Qoie^) ~ C'(6'). Consequently hairy black holes of charge 
~ 0{9) may be constructed in perturbation theory in 9 as in the previous paragraph (see Section 
[5]). The results of this analysis at small charge, combined with our numerical construction of 
these solutions at large charge, allow us to conclude that the lowest mass hairy black holes 
at fixed charge are never solitons. As in the previous paragraph, they appear to be possibly 
singular zero temperature finite entropy solutions. 

In summary, in this paper we have presented a rather complete picture of the solution space 
of a family of solitons, RN AdS black holes and a family of hairy black holes solutions. The 
solutions we have constructed certainly do not exhaust the set of static charged solutions in 
(1.1). For instance it was demonstrated in jl] that there exist an infinite number of 'excited' 
solitonic solutions (one based on each of the linearized excitations of the scalar field about global 
AdS^ space). Depending on the value of e^, some of these solitons may also be obtained as the 
lowest mass limit of a class of 'excited' hairy black hole solutions. At least at small charge, 
however, these new excited hairy black hole solutions are all unstable and so are unimportant 
for thermodynamical purposes. We feel that it is plausible that the thermodynamics of the 
system described by (1.1) is dominated by one of the three solutions studied in this paper 
at every value of mass, charge and e, although a verification of this suggestion would require 
further work. 

The results of this paper appear to throw up a numerical coincidence for which we do not 
have an explanation. The value = y appears to be special for two different reasons. First it 
marks the smallest value of the scalar charge for which one can construct arbitrarily small hairy 
black hole solutions. Second it appears, within numerical accuracy, to mark the dividing line 
for solitonic solutions. Below this value of e^, the solitonic branch of solutions constructed in 
this paper encounters the Chandrashekhar limit at large mass. Above this value of e^, smooth 
solitonic solutions exist at arbitrarily large mass. If these two different sharp changes do indeed 
occur at the same value of e^, this coincidence requires an explanation that we do not yet have. 
Of course it is entirely possible that the solitonic shift in Chandrashekhar limit does not occur 
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at = y but at a value that is coincidentally so near to this value that our numerics are 
unable to tell the two numbers apart. We leave a fuller consideration of this matter to future 
work. 

Another issue that we have not fully resolved in this paper is the precise nature of extremal 
hairy black holes. As we have reviewed above, the lowest mass hairy black holes appear to 
have zero temperature and finite entropy. Precisely at extremality these solutions appear to 
be singular (the curvature invariants and tidal forces grow large as we approach this configura- 
tion), however neither our perturbative techniques nor our numerical constructions are sensitive 
enough to determine what the nature of this singularity is. We leave a fuller characterization 
of these extremal solutions to future work. 

The nature of these extremal solutions also has bearing on the study of supersymmetric black 
holes in AdS^ x 5*5. As explained in f7j, IIB supergravity on AdS^ x admits a consistent 
truncation that sets equal the three diagonal U{l)s in the R symmetry 5*0(6), and includes 
a charged scalar field under this U{1). This consistent truncation admits hairy black hole 
solutions which were constructed perturbatively at small charge in |7]. At small charge these 
hairy black holes reduce to a supersymmetric soliton in the extremal limit. However the susy 
soliton does not exist beyond a critical charge of order unity. The nature of the extremal limit 
of hairy black holes beyond this charge is as yet unknown, and may have important bearing on 
our picture of the vacuum structure of A/" = 4 Yang Mills at finite 50(6) charge density. 

The plan of the paper is as follows. Section [2] provides a sharp summary of the phase 
diagram of the system (1.1) in the microcanonical ensemble (omitting some details discussed in 
later sections). Section [3j is entirely devoted to the study of the soliton family of solutions for all 
values of the scalar charge. Section |4] introduces the Reissner-Nordstrom— AdS (RN AdS) black 
hole; studies the Klein-Gordon equation for linearized charged scalar perturbations in the RN 
AdS background to find the onset unstable modes; and discusses the two sources of instability, 
namely the scalar condensation and the superradiant instabilities. Section |5] addresses the 
solutions with scalar field charge 3 < < y. It first uses a thermodynamic non- interacting 
model to find the leading order properties of the solutions and then it constructs these solutions, 
first using a perturbative approach and then a full numerical construction. Section [6] repeats 
this process but this time for solutions with scalar field charge > y. It starts with a 
thermodynamic non-interacting model to find the leading order properties of the solutions and 
then it constructs these solutions using a perturbative and a numerical construction. The 
details of the perturbative construction of the hairy black holes of section |5] and [6] are left to 
Appendices B, C and D. 



Note added: Reference 28 , which appears simultaneously with our work on the arXiv, 



discovers and studies a new branch of solitonic solutions in a model closely related to (1.1). 
The interplay of these new solitons with the solutions studied in this paper appears to be an 
interesting topic for a future study. We added a small section 3^ in a second version of our 
work to address some immediate questions. 



2 Summary of the phase diagram 



In this we section make the assumption that the thermodynamics of ( 1.1 ) is determined entirely 
by the three classes of static solutions considered in this paper. Under this assumption we 



present a summary of the phase diagram of the system (1.1) in the microcanonical ensemble. 
This phase diagram is qualitatively different depending on whether < 3, 3 < < y or 
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3 . We consider these three cases in turn. The phase diagrams presented in this system 
will be justified in detail in future sections. 
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Figure 1: Schematic phase diagram for < 3. Note that the sohton curve (green) is always 
above the extremal black hole (black). Pure RN AdS black holes exist for all values of masses 
above the extremality curve. The extremal black hole is never unstable and no hairy black hole 
solutions exist. The soliton exists up to a certain maximum charge, has a self-similar behaviour 
around a value Qcruie^) and is never the dominant phase. 



2.1 < 3 



At these values of e the phase diagram of (1.1) is always dominated by RN AdS black holes 



which are always stable. As the properties of RN AdS black holes do not depend of e, their 
phase diagram is independent of e^. At any given value of the charge the RN AdS phase exists 
down to a minimum mass; the phase boundary is given by extremal RN AdS black holes. 

Solitonic solutions exist up to a certain maximum charge, at these values of e, but they 
play no role in the thermodynamics of the system. In particular, solitons do not represent the 
ground state of the system at any given charge, as their mass is always greater than the mass 
of the extremal black hole at the same charge. 

The qualitative phase diagram for < 3 is shown in Fig. [Tj For a numerically constructed 
phase diagram with full details, see Fig. [sj (We find that an interesting cusp structure develops 
in the neighbourhood of Qcrit^ whose discussion we postpone.) 



2.2 3 < < f 

At these values of RN AdS black holes are unstable in a band about extremality for Q > 
Qo(e^). When they are stable, RN AdS black holes are always the dominant phase. New hairy 
black hole solutions are nucleated at the border of the RN AdS stability curve, and the system 
undergoes a second order phase transition to the hairy black hole phase on this line. Hairy 
black hole solutions exist down to masses below the RN AdS extremality bound, and represent 
the dominant phase whenever they exist. The lowest mass hairy black hole phase is extremal 
(zero temperature) and appears to have nonzero entropy. This solution is likely to be singular. 

As for < 3, solitonic solutions exist up to a certain maximum charge but play no role in 
the thermodynamics of the system. In particular, solitons never represent the ground state of 
the system at any given charge, as their mass is always greater than the mass of the extremal 
hairy black hole at the same charge. 
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Figure 2: Schematic phase diagram for 3 < < f . RN AdS black holes exist (red shaded) 
for all values of charge and mass above the extremality curve (black). Extremal black holes are 
unstable for Q > Qo(e^) and hairy black hole solutions exist for these values of charge. These 
solutions exist (blue shaded) between the curve of instability of RN AdS black holes (red) and 
a zero temperature hairy black hole solution (blue). Hairy black holes are the dominant phase 
whenever they exist. The soliton (green) exists upto a maximum charge Q = Qcru (where a 
cusp structure, to be discussed only later, appears) and is never the dominant phase. (The 
soliton curve can be below the extremal RN AdS line for Q > Qoie"^) but keeps above the 
extremal hairy black hole). 

Fig. |2]gives a qualitative picture of the phase diagram at these values of e^. For a numerically 
constructed phase diagram, see Fig. 



2.3 e2>f 

In this parameter regime, RN AdS black holes are unstable in a band about extremality, at 
every value of the charge. When they are stable, RN AdS black holes are always the dominant 
phase. As in the previous subsection, hairy black hole solutions are nucleated at the border of 
the RN AdS stability curve, and the system undergoes a second order phase transition to the 
hairy black hole phase at lower values of the mass. Hairy black hole solutions extend to masses 
below the RN AdS extremality bound, and represent the dominant phase whenever they exist. 

In this range of parameters, the lowest mass hairy black hole phase is an infinite temperature 
soliton for Q < Qczi^'^) (see Section [6] for details of the function (5c2(e^)). For Q > Qc^ie'^), on 
the other hand, the lowest mass hairy black hole is extremal (i.e. has zero temperature), has 
finite entropy, and is likely singular. 

The soliton exists at all values of the charge. As explained above, for Q < Qc^ the soliton 
represents the lowest mass hairy black hole or the ground state of the system. For Q > the 
soliton plays no thermodynamical role; in particular its mass is always larger than that of the 
lowest mass hairy black hole. 



Fig. 



3] gives a qualitative picture of the phase diagram for > y. For a numerically 



constructed phase diagram, see Fig. 24 
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Figure 3: Schematic phase diagram for > y. As for < y RN AdS black holes exist (red 
shaded) for all charges, and all masses above extremality (black). Hairy black holes exist (blue 
shaded) for all values of the charge, and the masses are below the curve of instability of RN 
AdS black holes (red). For Q < (5c2(6^), the lowest mass hairy black hole solution is a zero 



entropy soliton at infinite temperature (green) 
hole is extremal with finite entropy (blue). 



For Q > Qc^ie ), the lowest mass hairy black 



3 Solitons 

In this section we will numerically construct solitons at the full non-linear level. As we have men- 
tioned in the introduction, our system admits an infinite number of solitonic solution branches. 
The solitons we construct in this paper are those that are obtained by continuously increasing 
the scalar field at the origin, starting from AdS^ perturbed by the ground state linearized scalar 
fluctuation. We do not study any other branches of solitonic solutions in this paper. 

As we have explained in the Introduction, the space of solitonic solutions differs qualitatively 



depending on whether e > e^^^^^jj or e < e 



solcrit' 



As we have also remarked above, e 



solcrit 
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to within numerical accuracy. We first describe the methodology we use to construct our 
solitonic solutions, and then describe our results. 



3.1 Set up 

Following |4], we adopt the gauge: 

ds^ = -f{r)dt^ + g{r)dr^ + rUQl^^, A^dx^" = A{r) dt , = 0(r) , (3.1) 

where (/)(r) can be taken to be real wlog and dflj^^ denotes the standard metric on the round 
unit three-sphere. The equations of motion are (41: 



/2 



A" + 



f' + 



l + 2g( l + 2r 

2ge^r^ A^ 
3/ 



^2 ^ 
9 

2 1/2 

- r (p 



ge^ 
f 



0, 



A' -2ge^(f)^A = 



1 



l-g{l + 2r^)--r^(j)'^ 



f--{2ge'A'<l>'-A" 



(3.2a) 
(3.2b) 
(3.2c) 
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' 2 + — — + - r V" U + — 1 + 2 + — ^ = . (3.2d) 



r\ 3/3 / r \ 3/ 

where the prime ' denotes the derivative with respect to r. Notice that from the third equation 



in (3.2) we can ehminate g{r^ in terms of the other variables. Plugging this expression for g 



into the last equation in (3.2) yields a second order equation for /. We will take this equation 



together with the first and second equations in (3.2) as our set of fundamental equations to be 
solved. 

The boundary conditions that we shall impose are as follows. Since we want to consider 
asymptotically AdS^ solutions, the large r behavior of the various functions is given by 

A(r)=/i-^ + 0(l/O, 



e 



^4 



+ 0(l/r^) , (3.3) 



777 

/(r) = l + r2-- + 0(l/r^^ 



where [i is the chemical potential and e is the expectation value of the operator dual to the 
scalar field 0; the asymptotic charges of the solution are the mass M = and the charge 
Q = ^q. At the origin r = we shall impose smoothness, which implies a Neumann condition 
on all functions. 

For a given value of the scalar charge e, solitons form a one parameter family of solutions. 
This parameter can be taken to be the charge, for instance, but we find it useful to use /(O) 
instead. The reason is that, as discussed in [2], a necessary condition in order to have smooth 
horizonless solution is that /(O) > 0; in our gauge one has that /(O) = 1 corresponds to pure 
AdS and /(O) monotonically decreases as one moves along the soliton branch. Therefore, /(O) 
uniquely labels each solution for any value of e; on the other hand, we find that for e < Csoicrit 
the charge does not uniquely specify the soliton solution. 

The aforementioned boundary conditions can be easily implemented introducing a new 

2 

compact coordinate y = j^J^ and redefining the functions as follows: 

/(r) = (1 + r^) pf{y), 0(r) = e (l + r^) ~' p^y) , A{r) = p^iy) . (3.4) 

Then, the new functions j9j's should satisfy p/(l) = 1, Pa(1) = and p'^{l) = at infinity 
{y = 1), and at the origin {y = 0) we find: 

PfiO) = ^eVp,(0)V(0)^ pW = -"^^^^^^^ Pa{0) = \eh'pM'PA{0), 

(3.5) 

which follow from solving the equations of motion near the origin and imposing regularity. 



We have solved the equations (3.2) numerically using Newton's method and a Chebyshev 
pseudospectral collocation approximation. Alternatively we have also used shooting and the 
results of both methods agree. The data presented below was obtained using the pseudospectral 
method. 



3.2 Results: < e^i^^^, « f 

Our full nonlinear numerical construction of solitons agrees well with the perturbative construc- 
tion presented in fi] at small values of the charge (see for instance the right panel in Fig. [6]). 
The most striking qualitative feature of the solitonic branch of solutions, at these values of e^. 
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Figure 4: Left: qcru as a function of — e for < 32/3: note tliat tlie scale on ttie x axis 
is logaritlimic. The black dots correspond to our data. The red line is the fit to a linear behavior 
as in (3.6) in the region of small e. The blue line corresponds to a logarithmic fit of the form 

Qcrit = a In — + 6 of the data near e ~ \[^- This diagram suggest that qcrit blows up 

logarithmically as e approaches esoicHt- Right: Q as a function of /(O) for e = 3., 3.1, 3.2, 3.29, 3.5, 4. 
(from bottom to top in the plot). For e < esoicrit, Q tends to a finite value as /(O) — )• 0, and for 
e > esoicrit, Q diverges as /(O) 0. 



is that it terminates in a naked singularity at a critical value of the charge, (jcrit- This feature is 
readily explained at small e. Indeed, at very small e the soliton is simply an almost uncharged 
boson star. It seems intuitively reasonable that such a boson star hits the Chandrashekhar 
limit at a critical mass of order unity. At the first nontrivial order in e, the charge of the soliton 
may be measured by the gauge field A sourced by this almost uncharged boson star, and so is 
to order e. In other words to leading order in e, 

Qcrit = Ae, (3.6) 



Our full nonlinear numerical solutions, obtained using the methods of ^3.1 verify this expecta- 
tion, see Fig. |4] (left) and demonstrate that A ^ 0.0218. At larger values of e, Qcrit continues 
to grow. For e > e^^^^^j^ ~ y solitons exist at every value of the charge, i.e. Qcrit = oo. As we 

approach e = Csoicnt from below, the graph of Qcrit versus — ln(y^ — e) appears to asymptote 
to a straight line (Fig. Ill (left)) suggesting that the critical charge diverges like 




Qcrit ^ -0.0399 In W — - e + 0.159. 



We pause here to emphasize that we are inferring a divergence in Qcrit e = y y = 3.266, 
even though the largest value of Qcrit we have found in a simulation is 0.314 at e = 3.245. 
It turns out that simulations at values of e nearer to -i/ y are difficult; the slow growth of 
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Figure 5: Left: Kretschmann invariant evaluated at r = as a function of /(O) for the e = 1 
solitons. As /(O) — )• the Kretschmann invariant diverges at the origin, which signals the appearance 
of a curvature singularity there. The red dashed line corresponds to the value of the Kretschmann 
invariant for pure AdS^. Right: Phase diagram in the micro canonical ensemble for e = 1. On the 
y-axis of this plot we depict AM = M — Mext, where M^xt is the mass of the extremal RN AdS black 
hole with the same charge Q. RN AdS black holes occupy the shaded region and the soliton family 
of solutions is given by the black curve. This curve terminates at a naked singularity at some finite 
Q. In red we show the perturbative results of |;4!j; the agreement between the perturbative calculation 
and our numerical results is remarkable at small values of Q but they disagree at sufficiently large Q. 



the logarithm prevents us from obtaining more direct evidence for the divergence in tjcru as e 
approaches Csoicrit from belowj^ 

In order to further study the approach of e to esoicrit from below, in Fig. |4] (right) we have 
plotted Q as a function of /(O) for e = 3., 3.1, 3.2, 3.29, 3.5, 4. (from bottom to top curves in 
this plot). Recall that /(r) is the coefficient of —dt"^ in the metric of the solitonic solution; 
/(O) going to zero indicates a Chandrashekhar singularity in the solution. For e < esoicrit, /(O) 
vanishes at finite charge. As e approaches esoicrit the approach of the curves in Fig. |4] (right) 
presumably creep logarithmically up the y axis. For e > esoicrit the curves presumably never 
intersect the y axis. 



3.2.1 Approach to g^rit at e < esoicrit 

In the rest of this subsection we describe the approach of the solitonic branch of solutions to 
Qcrit at a fixed value of e < esoicrit in niore detail. 

At Q = qcrit the solution becomes singular. This singularity is signalled by the fact that the 



function / in (3.1) develops a zero at the origin r = at this particular (and finite) value of 
Qcrit- In turn, this implies that the curvature invariants diverge at this point. For instance, in 
Fig. Is] (left) we depict the Kretschmann invariant evaluated at r = as a function of /(O) for 



^In particular our data is consistent with the possibility that the linear behaviour in Fig. [i] levels out at very 
small e — esoicrit leading to a moderate finite value of (jmax at e = esoicrit- A definitive statement here needs 
further - preferably analytic - work. 
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Figure 6: q — Qcrit as a function of /o = /(O) for the e = 1 case. The black dots correspond to our data 



and the red curve is the fit given by the functional form (3.7). As we approach the singular solution, 
the regular solitons exhibit a (possibly) infinite set of damped self-similar oscillations similar to the 
behavior observed in j?]. 

the e = 1 case. Clearly, as /(O) — )■ this curvature invariant diverges and it can be shown that 
the divergence is like ~ /(0)~^. Ref. [i] had conjectured that the solitonic branch of solutions 
would cease to exist at a critical value of the charge; our results fully validate this conjecture 
in this range of parameters. 

As we approach the singular solution along the soliton family of solutions, the soliton ex- 
hibits a set of damped self-similar oscillations that (presumably) continue ad infinitum. More 
precisely, we find that near the singular solution the dependence of charge q of the solitons on 
/o = /(O) obeys a relation of the form, 

q - q,ru ~ a e-^'°^f^' cos (7 log f^^ + 6) , (3.7) 

for some constants (a, /3, 7, 5) that can be determined numerically from the data, and qcrit is 
the charge of the critical solution. We have analyzed different values of e (always with < y) 
and we have found that a and 5 do depend on e but /3 and 7 seem to be roughly independent 
of e. More specifically, we find (5 ~ 0.90 ± 0.05 and 7 ~ 1.20 ± 0.05 for the different values of 
e that we have checked. In Fig. |6]we illustrate this behavior for the e = 1 case. Therefore, for 
< y the phase diagram for the solitons (in the microcanonical ensemble) is the same as in 
the truncation of A^ = 8 gauge supergravity of [T]. We should emphasize that our analysis is 
fully numerical and it would be interesting to get some analytical results along the lines of [T]; 



in particular, it would be nice to calculate some of the constants in (3.7) analytically, get a 
better understanding of the singular solution and spell out the dependence on e (if any). 

In Fig. [7] (left) we depict the ADM charge Q vs. /(O) for the e = 1 case. We can observe 
that for small charges the perturbative results of (in red) agree very well with our numerics 
but they disagree at sufficiently large charges. In addition, this plot clearly shows that Q does 
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Figure 7: Left: Q as a function of /(O). Tlie cliarge exiiibits a (possibly infinite) series of damped 
oscillations around the critical value. Right e vs. Q. The expectation value of the scalar field as a 
function of the charge Q forms a spiral near the critical point. 



not uniquely parametrize the family of solutions. Instead, as discussed above, this quantity 
(and the other physical quantities too) exhibits a (possibly infinite) series of self-similar damped 
oscillations as we approach the singular solution, which corresponds to /(O) = 0. In the right 
panel of Fig. [7] we depict e vs. Q, which shows a spiralling behavior towards the singular 
solution. 

We close the discussion by noting that we have checked that the behaviour described in this 
subsection applies to values of that differ from y by less than 1%. We have also checked, 
on the other hand, that solitons exist for arbitrarily large charge for that exceed y by 1%. 
This is the basis of our claim that el^icru ~ y- 

3.2.2 Approach to the singularity 

In this brief subsection we study, from yet another angle, how the solitonic branch of solutions 
behaves as e is increased above e = Csoicrit- In Fig- [8] we display a set of graphs for the solitonic 
charge as a function of the value of the scalar field at zero, 0(0), for different values of that 
go through the phase transition. Note that for < e^^^^^jj the amplitude of the oscillations 
becomes larger as we approach e^^^^^jj, and the oscillations completely disappear for > e^^^^^j^. 
In other words, if we plotted 0(0) vs. Q we would see that the spiral unwraps for > e^^^^^j^. 
For > CgQj^^j^ we note that 0(0) approaches a constant for large Q and this constant is not 
too far from ^. This is consistent with a singular limit as in [29] and the results in Appendix [A| 
Similarly, the graph of Q vs. /(O) (see the right panel in Fig. [4| is monotonic for > e^^i^^^^. 

3.3 Results: > e^,^^^, « f 

For > y the behavior of the soliton family in the phase diagram changes completely. Whilst 
an unjustified extrapolation of the perturbative results of |4| to large charge suggest that solitons 
should have a Chandrasekhar bound for all values of e, our numerical results demonstrate that 
for > Csoicrit ~ ^ solitons can exist for arbitrarily large values of the charge. In Fig. M (left) 
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Figure 8: Charge of the sohton Q vs. (/)(0) for e = 1., 2., 3., 3.2, 3.29, 3.5, 4 (from bottom to top curves). 

"solcrit 



The behavior of Q as a function of 0(0) changes discontinuously as we vary across e^^i^^-^ w 32/3. 



The amphtude of the oscillations increases as we approach e'lgi^.^it from below and they completely 
disappear for > e^^^^^.^ . 

we plot /(O) as a function of the charge Q for the e = 4 case. The black curve corresponds 
to our data and the red curve corresponds to the perturbative results of |4|. Note that for 
small values of the charge the agreement between the perturbative results and our non-linear 
calculation is remarkable. On the other hand, an unjustified extrapolation of the perturbative 
results to large charge suggest |4j that /(O) should become zero at some finite value of Q, but 
our non-linear results show that this is not the case and in fact the soliton family of solutions 
seems to exist for arbitrarily large values of Q. The explanation for this is that for any > y 
the electrostatic repulsion is always strong enough to counter-balance the gravitational self- 
attraction, no matter how massive the soliton is. To check the regularity of the spacetime, we 
have computed the value of the Kretschmann invariant at the origin as a function of Q (see 
Fig. [9] right). Interestingly, for large values of the charge (and hence well in the full non-linear 
regime!) the value of the Kretschmann invariant at the origin reduces to that of pure AdS^ (red 
dashed line in the right panel of Fig. |9]), which implies that the curvature there is not large. 

As solitonic solutions exist at all values of the charge in this range of parameters, it is inter- 
esting to study the infinite charge limit of these solutions. In this limit the solitonic solutions 
should tend to solitonic black branes. Indeed we have numerically verified the following. Let 
L = g3. Then the limits 

gp{p) = lim L^g{pL) 

r . . y f{pL) 

fp[p) = hm —— 

Ap{p) = hm 
(j)p{p) = hm (j){pL) 

all appear to exist, implying that the large charge hairy black hole may be rewritten in new 
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e = 4 




Figure 9: Left: /(O) as a function of Q for tlie e = 4 case. In black we depict our numerical data 
and the red curve corresponds to the perturbative results of Q. The agreement between the 
two is excellent for small Q. An unjustified extrapolation of the perturbative results suggest, 
however, that /(O) should develops a zero at a finite Q. On the other hand our nonlinear results 
indicate that /(O) remains finite for any finite value of Q. Right: Kretschmann invariant at 
r = as a function of Q. This curvature invariant remains finite for any finite value of Q, and 
it approaches 40, (which is the value for pure AdS^) for large Q. 



coordinates r = Lp,t = T/L in the black brane form 

ds' = -fp{p) dr' + gp{p) dp' + p' {dx'f 
A = Ap{p)dT ^ ' 

In particular this implies that the mass and charge of our solitonic solutions scale, at large 
charge, like 

M^5Q^/^ M^Ce, (3.10) 



by conformal invariance. Indeed, as Fig. [TO] shows, our data approximately satisfies the afore- 
mentioned scaling behavior and yields B 0.9 and C ~ 1.0 at e = 4. The coefficients B and 
C are, of course, functions of e. 

While the solitonic solution is regular at every finite value of the charge, we have demon- 
strated in Appendix |A] that scaled solution (3.9) is not regular but instead has a singularity at 



r = at least at generic values of e^. This singularity appears, however, to be rather mild in 
nature. Neither the Kretschmann invariant nor (VRiem)^ of the scaled solution appear to blow 



up anywhere (see Fig. 11). It would be nice to better understand the near horizon behaviour 



of this planar solitonic solution along the lines of 29 ; however we leave this to future work. 

Later in this paper we study hairy black holes. As we demonstrate below, at finite but 
large values of the charge extremal black holes have lower mass than solitons. In the planar 
scaling limit described above, the solitonic mass approaches that of the hairy black brane, and 
it appears that the planar limit of the solitonic solution in fact coincides with the extremal 
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e = 4 



6 = 4 




logo 



Figure 10: Left: logM vs logQ for the e = 4 case. The black dots correspond to our data and 
the red curve corresponds to the best fit of our data to a function of the of form Af = BQ"' in the 



large charges regime. We find that a ~ 1.3, which is close to the scaling (3.10) required by conformal 
invariance, and B = 0.9. Right:. M vs. e. As this plot shows, for large values of these quantities 
the relation between them is approximately linear, in accordance with (3.10). We find M = Ce"^ 
C ~ a ~ 1.0. 



with 



hairy black brane. 
holes. 



We provide numerical evidence for this later, when we discuss hairy black 



3.4 



Second soliton branch for < e^^i^^-^ 



Reference 28 , which appears simultaneously with the first version of our work on the arXiv, 



discovers a new branch of solitonic solutions in a model closely related to (1.1). We confirmed 
with our numerical methods that this second branch exists for < el^icrit case, and that 

^"soicrity which agrees with the expectations 



it merges with the branch described above for 
28] Fig 



from 



mre 



12 



represents the two branches, for the particular value e 



3.2 < Ggolcrit- It 

is clear that one branch is continuously connected to small solitons, studied above, while the 
second branch exists for arbitrarily large charge/mass. We can study the large charge limit of 
the latter branch, as we did previously for solitons in the range > e 



In particular, the 

scaling behaviour (|3.10[), which follows from the expressions (|3.8[) and (|3.9[), should apply. We 



solcrit • 



confirm this expectation with our numerical analysis in Figure 13 



4 Instabilities of the Reissner-Nordstrom— AdS black hole 

We start this subsection with a review of the d = 5 Reissner-Nordstrom— AdS (RN AdS) black 
hole. We then study its instability with respect to the condensation of a minimally coupled 



Ref. [28j a lso finds that a second branch of solitons exists for > e^^i^^^^. Unlike the branch studied 



in section 3.3 which exists for all values of the charge, that second branch exists only for a finite range 

< Q min 

(e2) <Q<Q 



17 




Figure 11: Right: The Kretschmann invariant, K{p) for the scaled solution (3.9) may obtained from 
the Kretschmann invariant Ks{r) of the soliton as K{p) = Ymvi^oo Ks{r / L), where L = qs . We 
plot the RHS of this equation for 5 < g < 30. For the largest values of q shown in these plots the 
mass of the soliton is larger than 100 (in AdS units) and therefore these solitons are "large" and the 
considerations in ^3.3| should apply. Note that the curves fall nicely on top of each other (towards the 
right in the plot) and show no sign of a blow up at small p. Right: Note that the maximum value of 
the Kretschmann scalar asymptotes to a constant at large q, i.e., large Q. 

massless scalar, which signals the bifurcation to a hairy black hole. 

The RN AdS black hole is a two-parameter solution to the equations of motion that follow 



from (1.1), given by: 



V{r) ' V r V V 3 r 

/ 1 (4.1) 



0. 



where r = i? is the location of the event horizon and yU is the chemical potential. The existence 
of an event horizon requires that 

^'<lil + 2R'), (4.2) 
with the equality saturated by the extremal (zero temperature) black hole. 
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e = 3.2 



m 4 




Figure 12: Mass versus scalar field at the origin, witli two soliton branches represented for e = 3.2 < 



Ssoicrit- The upper branch, which exists for arbitrarily large charges, was found in 28 . The two 
branches merge for = e^^^^^^^. 



Here we record the basic thermodynamic quantities of the RN AdS black hole for later use: 



M 

Q 
s 

T 



RU1 + R' + 



3 R'^ 



2'= 2"^ 



TT 



(4.3) 



r;' 



2-nR 



2R' 



fx 



where M is the mass, Q is the charge, S is the entropy and T is the temperature of the black 
hole. The rescaled parameters m and q will be used later for convenience. 



4.1 Linear instability 



In this subsection, we consider the stability of the RN AdS black hole (4.1) under scalar 
condensation. Let us focus on the extremal solutions, 

V{r) =(^1-^^ (1 + + 2i?2) = 4 + _ i?)2 + o ((^ _ ^)3) ^ 

(4.4) 

IH / \ / / I \ 

At 



'Q{2+^]ir-R) + 0{ir-Ry) 



which are parameterized by the horizon radius R. The charge of the extremal black hole is 
given as a function of its radius by 



q = R\ -{1 + 2R^). 



(4.5) 
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Figure 13: Left: logM vs logQ for the second soliton branch, in the case e = 3.2. The black dots 
correspond to our numerical data and the red curve corresponds to the best fit of our data to a function 
of the of form M = AQ'^ in the large charges regime. We find that a ~ 1.287, which is close to the 
scaling (3.10) required by conformal invariance, and A ~ 1.095. Right:. M vs. e. As this plot shows. 



for large values of these quantities the relation between them is approximately linear, in accordance 
with (13. 101). The best fit to M = ^e" is a ^ 1.030 and A ~ 1.727. 



4.1.1 Near horizon instability for large black holes 



In the near horizon limit r — R <^ R, the metric (4.4) reduces to AdS2 x S^. In this region, the 
linearized equation for the charged scalar field (j) (see (1.1)) about this background reduces to 



the equation for a massive minimally coupled scalar with 

2 2 __3e^R^ l + 2i?2 

where lAdS2 is the radius of the AdS2 region. As is well known [l4] , a minimally coupled scalar 
in AdS2 is unstable whenever 

2;2 ^ 



''s''AdS2 



< 



It follows that the extremal black hole of horizon radius R is unstable whenever 

2(l + 3i?2)2 



3i?2(i + 2i?2)' 



(4.6) 



The RHS of (4.6) is a monotonically decreasing function of R. At large R, (4.6) reduces to 

1 



> 3 + 



2i?2 



+ C(l/i?4 



(4.7) 



It follows that very large extremal RN AdS black holes are unstable when > 3. The end 
point of the instability involves a condensate of the scalar field. By the Hawking area increase 
theorem it also has a horizon. Consequently, the end point of this instability is a hairy black 
hole. 
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It follows from the previous analysis that hairy black holes of charge q exist in the system 



(1.1) wheneve: 



(e (V9e2 - 24 - 3e) + 12) Je (3e + ^9e^ - 24) - 6 



?~ > ... . '..3/2 = ~^BF{e'). (4.8) 



24 (e2 - 3) 



This condition is not only sufficient, but also necessary for the instability of extremal black 



holes if -R > 1 (q' > ;^), as shown in |15j . 

4.1.2 Superradiant instability for small black holes 

For i? < 1, i.e for > y , the stability is not ensured by an analysis of the near- horizon region. 
Indeed, the leading instability of small extremal black holes is of the superradiant rather than 
Breitenlohner-Freedman varietj|^ (see [i]); the corresponding unstable modes are not localized 
entirely within the near horizon region. At small values of the black hole charge the extremal 
RN AdS black hole undergoes the superradiant instability provided 



?~>?~..(e^) = o\/o +^ 1-^ • (4-9) 




In summary we expect that hairy black holes exist in our system whenever 



q > go e , with S . , 2n / ~ / 2n r 2 ^ 32 ,\ ^'^^ 

go < q,r{e^) for > f {R <1), 



where the function go(e^) represents the smallest charge at which an extremal RN AdS black 



hole is unstablejj It was demonstrated in |4j that arbitrarily small extremal black holes in (1.1) 
suffer from the superradiant instability when > y . It follows that go(y ) = 0. Therefore it is 
natural to consider two different regions in the parameter space of hairy black holes, 3 < < y 
and > y. We will analyse these regions separately in later sections. 



4.2 Numerical results 

We now present a numerical analysis of the linearized instability of the RN AdS black hole, 
which complements the considerations in the previous subsection. We consider the wave equa- 
tion for a massless scalar field of charge e in the background of the RN AdS black hole and 
we look for marginally stable (that is, time-independent) spherically symmetric modes. The 
existence of such modes signals a bifurcation point in the phase diagram, and a new family of 
black holes with a non-zero scalar condensate (i.e. a hairy black hole) should emerge from the 
RN AdS family. We find that such modes exist provided that > 3, which indicates that hairy 
black holes will exist if that the scalar charge e satisfies this condition, as predicted above. 



solution to this equation into (4.5) 



In order to obtain this equation, we solved (4.6) to obtain R as a function of e , and then plugged the 



^This picture of superradiance can be reconciled with the more traditional version of an instability cj) ^ e~™* 
with onset mode w — efi hy a change of gauge. Consider the gauge A — A — fidt, for which the potential 
vanishes at infinity. A static scalar field in the gauge A (corresponding to the bifurcation to a hairy black hole) 
transforms to 4> ^ g-ie^tt ^j^g gauge A. 

*We assume here that all black holes of a given charge q are stable if the extremal black hole at that charge 
is stable. This expectation is intuitively reasonable, and has been borne out by all explicit computations to 
date. 
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The equation that we are going to solve is 

D^(p = (4.11) 

for 0(r) real and where = — i e is the gauge covariant derivative on the RN AdS back- 
ground. As discussed above, we will only consider time-independent and spherically symmetric 
modes and, furthermore, we will impose the following asymptotic behavior, 

0(r) ~ for r-j-oo, (4.12) 

so that e is the vacuum expectation value of the operator dual to 0. In addition, we will also 
require that is regular at the horizon. These boundary conditions can be easily implemented 
redefining the field 0(r) as 

^{r)=p4r)(-] , (4.13) 



with 

p;Loo = 0' (p;-4p,)U = 0. (4.14) 
With these boundary conditions we can solve ( |4.11[ ) as a boundary value problem using shooting 



as in 8 . However, in this paper we will follow a slightly different route. Following 15 



we cast 



(4.11) as a generalised eigenvalue problem, 

L(r)0(r) = A(r) 0(r) , (4.15) 

where the scalar charge e appears as the generalised eigenvalue. Here L{r) is a second order 
linear differential operator. Then, for a given RN AdS background uniquely specified by {R, fi), 
our strategy consists in finding the eigenvalue e for which there exists a mode that satisfies the 
above boundary conditions. 



The results are depicted in Fig. [Mj For a given black hole size R, we find that the minimum 
value of for instability is obtained for black holes in the extremal limit. Also, our results 
indicate that the minimum value of monotonically decreases from y to 3 as the black hole 
size increases, and for < 3 all black holes are stable under scalar condensation. On the other 
hand, for > y all extremal black holes are unstable. 



These results agree with the analytical predictions. In Fig. 15, we plot the minimum charge 
of unstable RN AdS black holes for values of close to (but below) y; the prediction of (4.10) 
is confirmed. 



5 Hairy black holes for 3 < < ^ 



When < 3, solitons and RN AdS black holes are the only static charged solutions of the 
Lagrangian (1.1). When > 3, however, the spectrum of charged static solutions also includes 



hairy black holes as the linearized analysis in ^4.1 and ^4.2 suggests. 

In this section we will investigate the spectrum of hairy black hole solutions in the parameter 
range 

3<e2<— . (5.1) 



32 

y 



While our investigations are partly numerical, we have also been able to obtain some analytic 



results at the upper end of the parameter range (5.1) as we now explain. 
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Figure 14: Critical value of the scalar charge e for the existence of an instability as a function of the 
mass M and charge Q of the RN AdS background. The minimum value of monotonically decreases 
from ^ to 3 as the size of the black hole increases. AM is the mass difference with respect to the 
extremal RN AdS black hole of the same charge. 



Whenever (1.1) hosts hairy black hole solutions of parametrically small charge, these so- 
lutions may be constructed analytically using perturbative techniques of Ej. As go(e^) is of 



order 6 when = y (1 — 6'), it is possible to use the perturbative techniques described above 
to completely characterize the spectrum of hairy black holes with charges of order 9 at these 
values of e^. In the first two subsections below we present an analytic construction of hairy 
black holes, perturbatively in 9. We use the results of our analysis to conjecture the qualitative 
properties of hairy black holes in the full range 3 < < y. In the rest of this subsection we 
present the results of a detailed numerical analysis that give evidence for our conjectures and 
fill out several details. 



5.1 Non- interacting model for small hairy black holes at small 9 

In this section and the next, we study small hairy black holes for values of scalar charge 
close to but less than y. We set 
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e' = ^{l-9) (5.2) 



where ^ is a small positive number. In ^5.2 we demonstrate that our system admits a spectrum 
of hairy black hole solutions at small 9 by explicitly constructing the hairy black holes in 
a perturbative expansion in 9. This construction is rather involved but the final result for 
the thermodynamics is rather simple at leading order in 9. It turns out that the leading 
order thermodynamics of these small hairy black holes is reproduced by modelling them as 
a non interacting mix of a RN AdS black hole and a soliton (this is in accord with previous 
experience in |4] and [?]). As a prelude of, and motivation for our perturbative construction, 
in this subsection we work out this simple thermodynamical model of small hairy black holes. 
The formulas obtained in this subsection will be derived more systematically (as the first terms 
in a series expansion) in the next subsection. 
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5.1.1 Thermodynamics of small RN AdS black holes 

Small RN AdS black holes appear in a two parameter family labeled by their radius R and 
their chemical potential ^bh- In this paper we will find it convenient to parameterize R and 
liBH in terms of auxiliary variables a and a defined by 

R^ = ae. 



2 

f^BH 



[1 + aaO) 



(5.3) 



where a and a are positive numbers. In this subsection and the next we are interested in the 
small 9 limit with a and a of order unity, a parametrizes the radius R of the black hole in units 
of 6 and a measures the deviation, in units of 6', of the chemical potential squared of the black 
hole from a critical value |. The usual formulas of black hole thermodynamics determine all 
thermodynamical properties of the black hole as a function of a and a; 
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^^BH — 


- (1 + aaO) 
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niBH = R^ 



1 + R' + -iiIh 



Qbh 
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(5.4) 
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where 



(5.5) 
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Note that all nonsingular black holes have a < 2. Black holes with a = 2 are extremal. Black 
holes with a > 2 formally have negative temperature and are unphysical (they have naked 
singularities). Note also that when a and a are of order unity, the black hole mass and charge 
is of order 9. 



5.1.2 Thermodynamics of small charge solitons 

As we have explained above, in addition to black holes, our system admits regular solitonic 
solutions at all values of e^. In [i] these solutions were constructed in a perturbative expansion 
in their charge. This construction allowed a determination of the energy rrisoi and the chemical 
potential figoi of the these solitons as a function of their charge qsoi at small charge. In this 
section we are interested in solitons whose charge is of order 6, and so we set Qsoi = bO where 
the small number 6 was defined in (5.2) and 6 is a positive number of order unityj^ Plugging 
into the formulas of [4] we find 

Qsoi = b9 , 

f^soi = l+(j-^^ qsoi + O(gL) =2qc + e (^ge + y ) + 0(9') , 

16 2 / 64\ 2 3 ^ (5-6) 

\2q, 7 J V 7 339570qJ ^ ' 

As the solitonic solutions have no event horizon they have zero entropy. As these solutions can 
be continued to Euclidean space with a thermal circle of arbitrary size, their temperature is 
indefinite. 



5.1.3 Hairy black holes as a non-interacting mix 

Consider a system with net charge 

q = q9, (5.7) 

and net mass 

m= -9+\-q^ + Am ) 9'^ . (5.8) 

The parameters q and Am parameterize the hairy black hole we wish to study. In this section 
we assume that q is held fixed as 9 is taken to zero. It will turn out that hairy black holes 
with charge q9 occur over values of Am that range from approximately to a negative number 
independent of 9. In other words Am is always order unity or smaller. 

In this subsection we will model hairy black holes as a non interacting mix of solitons and 
RN AdS black holes. A hairy black hole can partition its charge and mass arbitrarily between 
the black hole and soliton phases, and chooses to do so in the manner that maximises its 
entropy. 

How does a system maximise its entropy? As the soliton carries no entropy, all of the 
entropy lies in its black hole component. Let us suppose that the total mass m is partitioned 
up between black hole and soliton as m = m^o/ + TnBH and a similar partitioning for its total 
charge q = qsoi + qBn) 

St = SBniqBH, rriBH) + Ssoi{qsoi, rrisoi) = Ssniq - Qsoi, m - rrisoi) ■ 
^Note that while the small number is a parameter of the theory, b parametrizes the solution we study. 
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Maximising St w.r.t. qsoi and use of the first law gives 



dm 



sol 



dSBH 
dlUBH 



dq. 



sol 



IJ-BH — l^sol 



(5.9) 



In other words we must partition up charges so as to ensure that the chemical potential of the 
black hole equals that of the soliton. Assuming the model spelt out above, we will now describe 
the spectrum of hairy black holes in our system. We will first work this out at leading order 
in It turns out that the leading order result has a degeneracy that is lifted at next order in 
6] for this reason we will also compute the next correction in 6' to a particular aspect of the 
thermodynamics of the non interacting model. 



5.1.4 Hairy black holes at leading order in 9 

Let us suppose that the system divides itself into an RN AdS black hole parameterized by a 
and a and a soliton parameterized by b. We will now determine a, b and a in terms of q and 



Am (as defined in (5.7) and (5.8)). To the order of interest in this subsubsection 



ruBH = ^aO + a^(l + a)e^ , 
qBH = '^aq^e + aa^q^'^ + 0{f) , 



(5.10) 



while for the soliton one has 



fJ^soi = 2qc + 9 {qc + 



msoi 



-9 



2g. 



+ 
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y 



Qsoi = bO . 

Charge and mass conservation imply that 



q = qe = qsoi + qsH ^ q = i^aq^ + b) + 0(6) , 
q. f2 



7 



(5.11) 



m 



-9 +{-q' + Am] 9' 
qc \3 



rrisoi + rriBH Am 



b 

2^ 



ab 862 
^ ~ 21 



0{9). 



Eq. (5.12) gives 



2qc 



+ 0{9) . 



(5.12) 
(5.13) 

(5.14) 



and (5.11)) 



The requirement of chemical potential matching (5.9) gives (using the first equations of (5.10) 

36 



36 

y 



+ 0(9) 



(5.15) 



where we have used (5.14) in the last step. Plugging (5.14) in (5.13) gives an equation for 6. 
Solving this, we find 



7 
10 



(g - qc) ± 



7 

To' 



1 2 



21 A 

H Am 

20 



0{9). 



(5.16) 
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Recall that the black hole component of the mix must have a < 2. From (5.15), however, this 



implies that b < Tr{{Q — Qc) + C){6). As 6 > 0, it follows immediately that no solution exists for 



q < qc + 0{6). It also follows that acceptable roots for b (to leading order)in (5.16) are given 
by 



7 

To 



{q-qc 




+ 0{e) when q>qc + 0{e). (5.17) 



This solution is physical only when 6 > 0; this requires 

Am < 0{e). 

Moreover b must be real; this requires 

Am>-^iq-q,f + Oie). 

Note that a = 2 when the last inequality is saturated. It follows that the black hole component 
of a hairy black hole is extremal at the lowest allowed value of Am at any q. 
In summary, we have a solution whenever 



and 



q>qc + 0{e) 



7 



0{e)>Am>--{q-q,y + 0{e). 
15 



(5.18) 
(5.19) 



This range is plotted in Fig. [T6| In this range the black hole component of the mix has 



a 



2qc I 10 10 



1 / 3g 7 _ 3 
a I 10^ ^ 10 ~ Tg^; 




(5.20) 



10 



(q-qc 



21 

+ -Am I + 0(6) . 



The entropy, temperature and chemical potential of the hairy black hole are given by the 



formulas (5.4) with the values of a and a in (5.20). Note that 



- 2qc\l0 10 J ^ ' - 2 



(5.21) 



where we have used (5.18) in the last inequality. As we have explained above, at a fixed value 



of q, Am varies in the range (5.19). As Am is lowered within this range the hairy black hole 



horizon area - parameterized by a - decreases while its chemical potential - parameterized by a 
increases. The temperature of the hairy black hole also decreases upon lowering Am, reaching 
zero at the lowest allowed value of Am. 



5.1.5 Splitting the degeneracy at 0{9) 



The phase diagram depicted in Fig. 16 has the property that RN AdS black holes and hairy 
black holes coexist only on a single line, namely Am = 0. This odd feature is an artefact of 
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Figure 16: Range of allowed values of Am for hairy black holes (blue shaded). Note that hairy 
solutions exist only for q > q^. The upper limit, Am = denotes the onset of superradiant 
instabilities and corresponds to 6 = (red line). This is the also the extremality line for the 
pure RN AdS black hole. The lower limit, Am = —j^{q — qcY corresponds to the extremal hairy 
black hole with a = 2 (blue line). RN AdS black holes exist for all values of q and Am > 
(pink shaded). The soliton hes on Am^oZ =9^ (green curve). 



working at leading order in the 6 expansion. The actual phase diagram includes a region of 
coexistence of these two phases. It turns out that the height of this region is of 0{9) (in the 
variable Am) as we will now explain. 

The extremality curve for extremal RN AdS black holes is given by (see |4J) 



Am. 



ext 



3 V 3 



(5.22) 



We will now show that hairy black holes start existing at larger values of Am, demonstrating 
that hairy black holes and RN AdS black holes coexist over a range of charges and masses. The 
instability curve for RN AdS black holes is given by equating the chemical potential of these 
black holes to f (= 2qc + qcO + 0{e'^))^ This is the condition (from (|K9|) and (Q) 



a 



+ o{e). 



(5.23) 



Plugging (5.23) into the second and third equations of (5.4), we get 

q =2aqc + 0(6) . 



Eliminating a from (5.24) and (5.25), we get 

2 

Amsr = 



Qc 
2 



e + o{e'^ 



(5.24) 
(5.25) 

(5.26) 
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-1.5L 



I g^^ (yellow dashed line) 



Figure 17: Plot of Anigr = -y - f ) d{red line), Amext = - 
and the allowed range Am for hairy black holes (blue shaded) for 6 = 0.03. Hairy black holes 
exist for all q > Qc and Am between the red curve (onset of superradiance) and the blue curve 
(extremal hairy black hole). The yellow dashed curve denotes the extremality line for RN AdS 
black holes. Pure black hole solutions exist above this line for all values of q. RN AdS and 
hairy black holes coexist in the region between the solid red curve and dashed yellow curve. 



A plot of Amsr and Amext is shown in Fig 17 Note that for q < qc, Amsr > Amext- In other 
words, the RN AdS black holes reach their extremal limit before the onset of superradiant 
instabilities and therefore, hairy black holes do not exist for these values of charges. 

For q > qc, Amsr > Amea;|^ For these values of charge q, the RN AdS black holes suffer 
from a superradiant instability before (i.e. at a higher mass than) they go extremal; the region 
of the phase diagram between these two masses therefore hosts both RN AdS and hairy black 
holes. 



5.2 Perturbative analysis of hairy black holes 



In §5.1 , we have presented a thermodynamic model that makes the following predictions for 
small hairy black holes in the theory with 
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1. Hairy black holes at charge q exist only when 

q>qc9 + 0{e^), 



2. We study hairy black holes with charge q ~ 0{6), and so define q 



(so that q is of 



^"This is the minimum value of chemical potential of the soliton. 
^^The inequality is saturated at q = qc- 
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unit order at small 6) . We also parameterize the mass m of hairy black holes by 

m = 9- + 9^-q^ + Am] . (5.27) 
The thermodynamic model predicts that hairy black holes exist only in the range 

(g - I) + 0(9') >Am>-^{q- q^f + 0{9). 
Note in particular that Am is always of order unity or smaller at small 9. 

3. The lowest mass hairy black hole is roughly approximated as a non-interacting mix of 
an extremal RN AdS black hole and a soliton. It has zero temperature, but non-zero 
entropy. 

4. The value of the scalar field in this solution is is of 0{\/9). Hairy black holes are approx- 
imately a non interacting mix of a RN AdS black hole of radius R ~ a/^, charge q ^ 9 
and energy above extremality ~ 9^. 

In this section we justify all these predictions by presenting a perturbative construction of 
hairy black holes at small 9 and charges of order 9. 



5.2.1 The mathematical problem 

In this subsection we describe the mathematical problem we require to solve in order to con- 
struct hairy black holes in our model. 



The equations of motion for the Lagrangian (1.1) are 



0, 



rpEM _i_ rjimat 



where 



(5.28) 
(5.29) 

(5.30) 



T. 



EM 



TJZ"' = {D,<P){D,<py + (D^0)*(D.0) - g,^D,^ 



We are interested in stationary, spherically symmetric solutions of the system described by ( 1.1 ). 



As explained in ^3.1 we adopt a Schwarzschild-like gauge (3.1). The four unknown functions 
/(r), (?(r), A[r) and 0(r) are constrained by Einstein's equations, the Maxwell equation and 
the minimally coupled scalar equation. With this ansatz, the equations of motion boil down 
to the equations (3.2). It was shown in [4] that these equations admit a 6 parameter set of 
solutions, 
A{r) = 4 



One of the solutions is empty AdS^ space, given by /(r) = + 1, g{r) = 
(r) = 0. We are interested in those solutions that asymptote to AdS spacetime, i.e. 



solutions whose large r behavior is as in eq. (3.3) 



/(r) =r^ + l + 0{l/r^) 

g{r) = -^,+0{l/r'] 
1 + 

A{r) = 0{1) + 0{l/r^) , 
0(r) = C(l/r^) . 



(5.31) 
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It turns out imposing these conditions on / and eliminates 2 of the 6 solutions so that the 
system of equations admits a four parameter set of asymptotically AdS solutions. We will also 
be interested in solutions that are regular in the interior. In particular we demand that the 
functions / and do not blow up anywhere apart from at the singularity of the solution (so 
that our solution has a smooth event horizon). These two requirements generically cut down 
solution space to distinct classes of two parameter solutions In this section we will determine 
the two parameter family of hairy black hole solutions at small 9. 



5.2.2 Parametrization of solutions 

In the next subsection we will describe a perturbative construction of hairy black hole solutions 
at small 9. In this subsection we will specify how the perturbative solutions we construct are 
parameterized. We use the parameters k and a to label solutions, as we now explain. 
k is defined by the requirement that 



h9 

+ 0{1/t') (5.32) 



^4 



at large r. 

a parametrizes the radius of our hairy black holes. Let R denote the radial location (in the 



coordinate system of (3.1)) of the event horizon of the hairy black hok^^ The parameter a is 
then defined by the relation 

R^ = a9. (5.33) 

We wish to construct hairy black holes with k and a taken to be of order unity, in a 
perturbative expansion in 9. As is clear from the previous section (and as we shall see below), 
the solutions we will construct are approximately a non interacting mix of an RN AdS black 
hole and a soliton. The RN AdS component of this mix has radius of order \/9, charge and 
mass of order 9. The soliton component of this mix has charge and mass both of order 9. 

The hairy black hole is constructed in a perturbation expansion in \/9. The starting point of 



this expansion is an RN AdS black hole labeled by a and a (see subsection 5.1.1 ). This solution 
is then perturbed by a scalar condensate in order to meet the requirement (5.32). It turns 
out that the three parameters k, a and a cannot be independently varied. The requirement 
of regularity of our solutions fixes a as a function of a and k. To leading order, the relation 



so determined is in precise agreement with the results of the non interacting model of ^5.1 
justifying the non interacting model. 



5.2.3 Nature of the perturbative expansion 

In this subsection we briefly explain the nature of the perturbation theory we employ in order 
to construct hairy black hole solutions in a perturbative expansion in 9. Many more details of 
our procedure, together with a detailed listing of all our results are presented in Appendix |B] 

\/9 parametrizes the amplitude of the scalar field perturbation about the RN AdS black 
hole that constitutes the starting point of our perturbative expansion. However 9 also appears 
as a parameter that governs the size of the starting RN AdS solution. 

In the 'far field region' r ^ R = \/a9 where the RN AdS solution is a small perturbation 
about AdS space, the expansion of our solution in 6* is a standard expansion in the amplitude 
of the scalar perturbation about unperturbed AdS space. When r ~ \/9 the starting solution 
of perturbation theory, i.e. the RN AdS black hole, depends on ^ in a crucial way. However it 

^^More invariantly R is defined so that the area of the event horizon is given by Ah = . 
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turns out that the dependence of the RN AdS solution on 6 in this region is rather simple, and 
may be scaled out of a problem with an appropriate choice of coordinates and variables. As a 
consequence, in this region as well, the expansion of the hairy black hole solution in 6 is simply 
an expansion in the amplitude of the scalar perturbation, albeit in the appropriate coordinates. 

In summary, the perturbative expansion of the hairy black hole solution in 6 is everywhere 
an expansion in the amplitude of the scalar field, though we need to employ different variables 
in different regions in order to make this apparent and useful. In order to practically implement 
our perturbative expansion, we found it useful to work separately in three regions. In the far 
field region r 3> -R we worked with the usual AdS radial variable r. In order to describe the 
other regions of interest to the perturbative analysis, we need to introduce some new notations. 
Let Rin denote the inner horizon radius of our RN AdS black hole, while R is its outer horizon 
radius. In the parameter regime of interest to this paper it turns out that R—Rm ~ \CR^ ~ ^^^^ 
i.e. the black holes we study are very near extremality) . The intermediate field region is 
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defined by the condition r -C 1 but r — i? 3> CR^- this region our solution admits an 
expansion in 9 when expressed as a function of the scaled variable y = j^. In the so called near 
field region, r — R <^ R — Rin, on the other hand, our solution admits an amplitude expansion 
in 6 when expressed as a function of the variable z = = In the small 9 limit 

of interest to this subsection, the far field and intermediate field regions share a substantial 
overlap. Similarly the intermediate and near field regions overlap substantially. For this reason 
we are able to construct the full solution by working individually in the three region and then 
patching the full solution up using a matching procedure (following [i] and [7|). 

In Appendix |B] we have discussed the general structure of our perturbative expansion in 
some detail. In subsection |B.l of that Appendix we have described in detail the general method 



employed in our perturbative expansion. In subsections B.2 and B.3 we have described the 
explicit implementation of this procedure at order \/9 and 9 respectively. Finally in subsection 
we present a detailed listings of the final solutions of our hairy black holes to 0(9^^"^) in 



B.4 



each of the far field, intermediate field and near field regions. In the rest of this subsection we 



compute the thermodynamics of the solution presented in subsection B.4 



5.2.4 Thermodynamics 

It is a straightforward exercise to compute the thermodynamical charges and potentials of the 



black holes presented in Appendix B.4 For convenience we work with the rescaled mass and 



charge variables q and m defined and employed in j4j 

Q = 

i (5.34) 
M = — m. 
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where C = 2 — ap is independent of 9. 
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The formulas we have employed to compute the mass and charge and other thermodynamic 
quantities from the solutions of Appendix B.4 are 



T 
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I out I 



/i = lim 

r— >oo 

Ti= lim [r^ (l + r2-/°"*(r))l 

1 — s-oo 

q= lim [r^ (/x - A°"*(r))] , 



47rC (a^)^/^ 



r{z) 



(5.35) 



2=0 



(see Appendix B.l for a definition of the functions /°"*, and /*"). Employing these formulae 
we findo 
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(5.36) 



We have verified that these quantities obey the first law of thermodynamics 

dM = TdS + fidQ . 



^Throughout this paper, we follow ^ and consistently omit the factor of G5 ^ from all our extensive quantities 
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5.2.5 Comparison with the thermodynamics of the non-interacting model 



In this subsection, we demonstrate that the thermodynamical formulae (5.36) reduce to (5.12) 
and (5.13) at leading order. Using (5.36) and definitions (5.7) and (5.8), we can calculate q 
and Am for this system. Upto the order of interest in this section, 

(9a + Ak) 



Am 
T 



3^6 
Ak 



+ oie) 



567 
(42a - 8A; 



126a + 16A; - 63) + C(^) 
211 



(5.37) 



A27ry/a 



In order to make contact with the results of §5.1[ we need to identify the parameter k defined 
in (5.32) with the parameter b of ^5.1, We do this as follows. 
The pure soliton with 



(r) 
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has charge 



Qsol 



3q, 



(see Q). From the definition of b (see (5.6)), we immediately conclude b 
terms of a and b the formulae ( 5.37[ ) are modified to 



k 
3qc 



Am 



q = {2aq, + b) + 0(6) 
b ab 

2^"^~ 2r 



+ 0{9) 



Written in 



(5.38) 



These are exactly the formulae (5.12) and (5.13) of the non-interacting model. 

At leading order, the hairy black hole can therefore, be thought of as a non-interacting mix 
of an RN AdS black hole and a soliton. 

We now compute the equation for the line denoting the onset of the superradiant instability 
at one higher order in 6. This line is obtained by setting = in (5.36). To the order of 
interest (in terms of q and Am) we have 



Am 



q = 2aqc + 0(6^ 
1 



a a9 + Oie'). 



(5.39) 



Eliminating a from (5.39), we get 



Am = q (q 

3 



2 



9 + 0{9^ 



(5.40) 



as predicted by (5.26) in ^5.1.5 
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5.2.6 Break down of perturbation theory very near extremality 

Note that our spectrum of hairy black holes includes a one parameter set of extremal hairy 
black holes; at leading order this occurs when 42a - S/c - 21 = 21a C = (see ( |536| )) and 



in particular when ( = 0. While all thermodynamical formulae in (5.36) are smooth in this 
extremal limit, the solutions themselves develop a singularity, as we now explain. 

The small hairy black hole solutions we have constructed in Appendix |B] are regular (away 
from the black hole singularity) . In particular the value of the scalar field is finite everywhere 
outside the outer event horizon; the value of the scalar field at the event horizon is of order \/6 



Strictly at extremality, however, the scalar field diverges at the event horizon (see (B.72)). 
As our perturbative expansion is essentially an expansion in the amplitude of the scalar field, 
this divergence indicates a breakdown of our perturbative techniques (and so in particular is 
not a self-consistent and reliable prediction of our perturbation theory). It really indicates only 
that our perturbative construction of hairy black holes does not apply all the way down to 
extremality. 

How near to extremaliity does perturbation theory break down? We can examine the 



breakdown in perturbation theory in more detail as follows. As is clear from (5.36), all hairy 
black holes constructed in this paper have T ~ y/O. The parameter t = ^ characterizes the 
departure of our hairy black holes from extremality. 

^ _ T _ y/^C 



It turns out (see (B.72)), that at very small values of t the value of the scalar field at the horizon 

e expans 

^3/2 Int. 



is proportional (at 0{9^^'^) in the perturbative expansion) to 

4a 



k 

This term becomes comparable to the scalar field at leading order (signalling a breakdown of 
the perturbative expansion) when 

k (5.41) 
Consequently perturbation theory breaks down at exponentially low values of the temperature. 



Our perturbative construction of hairy black holes, and so the thermodynamical formulae (5.36 ) 



and the results of the non interacting model, all strictly apply only above the temperature 



(5.41). In particular our perturbative results cannot reliably be used to study extremal hairy 



black holes, although we can get exponentially near to these solutions at small 6. 

5.3 Numerical results 
5.3.1 Method 

In this subsection we construct the hairy black holes at the full non-linear level by numerically 
solving the equations of motion. Our ansatz for the metric, gauge field and scalar field are 



the same as in Eq.(3.1), and therefore the equations of motion are given by (3.2). We shall 



require that the asymptotic behaviour of the various fields is the same as in the soliton case. 



i^C = 2 - ao = (see footnote 13). 
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which is given in Eq.(3.3). Of course, the main difference now is that we have to require that 
the spacetime has a (regular and non-degenerat^Q event horizon; this imphes that /(r) has a 
simple zero at r = i?. The parameter R then determines the radius of the horizon. Furthermore, 
we choose a regular gauge for the gauge field A on the horizon, such that At{r = R) = 0. 

To implement the boundary conditions we find it useful to introduce a compact radial 
coordinate, y = 1 — \ and redefine the functions /, and A as follows. 



with r 



A{r) 



/=^. The boundary conditions that we shall impose are: 



{1 + r')-' p^{y) 

^ PA{y) 



(5.42) 



at infinity and 

p'/(0) = 



12 



(8 + e2p^(l)2) =0^ 



(5.43) 



p;(o) 

Pa(0) 



- 6j9^(0) (1 + R^) (1 + R^f {3pf{0)R'' - 2pA{0y (1 + 3i?2)) 

+ e\%{0)^pA{0)^R^ {2pa{0)^ + SpfiO) (l + R^)) 

X "l8 (1 + R^Y (pfiO) + 2pjiO)R^) (1 + R^ 
2p0(O)i?4 



-1 



(5.44) 



e'e'p^{0)'pA{0)R' {2pA{0f + 3pf{0) (1 + R')f 
36pf{0y (1 + 2i?2) (1 + + ^4 + ^6)2 



The latter boundary conditions follow from simply solving the equations of motion in a near- 
horizon expansion. 

We construct, in the present range 3 < < y, hairy black holes whose radius is not too 
small compared to the radius of AdS^, but we find no evidence for the existence of very small 
hairy black holes. This is in agreement with the previous perturbative analysis. Instead we 
find that the extremal limit of the hairy black holes is singular. It does not coincide with the 
soliton branch studied, as we also predicted from the perturbative treatment. 



5.3.2 Results 



The results are presented in Fig. 18, On the left, we have plotted AM vs. Q for e = 3.2 
{6 = 0.04) for both the soliton and the hairy black holes; AM is the mass difference with respect 
to the extremal RN AdS black hole of the same charge. We have added the analytical predictions 



from the perturbation theory, based on (5.27) and (5.36), and there is good agreement for small 



charge as it should. On the right, we have also plotted the same data with the entropy as third 
axis. Other values of e within the range 3 < < y give a qualitatively similar picture. As this 



16 



As mentioned in the Introduction, regular extremal hairy black holes are not allowed in the system under 



study 27 . 



36 



Fig. shows, there are no hairy black holes with arbitrarily small charge. As an extra check of 



our results, we represent in Fig. 19 hairy black holes along the merger line (upper mass bound) 
and along a line of constant boundary condensate e = 0.1, comparing the numerical data with 



the prediction of (5.36) for small charge. 



It is clear that the soliton family does not arise as a zero size limit of the hairy black hole. To 
investigate this, we plotted the temperature T and the Kretschmann invariant evaluated at the 
horizon, as a function of the horizon size, keeping e = 0.5 fixed; see Fig. |20j The plots suggest 
that the minimum size of hairy black holes corresponds to an extremal (i.e. zero temperature) 
limit which is singular. Since very close to extremality our numerics are not reliable, we 
have used extrapolation to estimate the value of R for which the black hole would have zero 
temperature, and we find that T = for R ^ 0.162. Similarly we find that the Kretschmann 
invariant would blow up for R 0.154; which is close enough, given the extrapolation. Finally, 
as predicted before, notice that the solitons are more massive than the extremal hairy black 
holes of the same charge. 

For other values of the scalar condensate e in this intermediate e range, the picture is 
qualitatively the same. 
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Figure 18: Left: AM vs. Q for e = 3.2. The shaded region is the area occupied by RN AdS black 
holes, the black curve corresponds to a soliton branch and the blue region corresponds to hairy black 
holes (the 'horizontal' lines of the blue grid correspond to fixed values of e, and R decreases to the 
left; the 'diagonal' blue lines near the green curve are segments of a hairy black hole with fixed R, 
and e grows to the right). The lower mass bound of hairy black holes is well described at small charge 
by the dashed green line, which is the perturbative prediction. The red curve is the line of marginal 



modes of the linear problem discussed in section 4.2 it agrees with the dashed magenta line for small 



charge, which is the perturbative prediction. It is clear that the soliton curve and the hairy black 
holes surface are not related in the range 3 < < 32/3. Right: same data with the entropy as third 
axis. The red surface is the strip of RN AdS black holes between the marginal line of stability and 
extremality. 



In this section we have described the construction of a hairy black hole solutions labelled, 

4 n 

for instance, by their mass Q and charge M. Let us now set M = (Q^. If Q = L is taken to 
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Figure 19: These plots M vs. Q represent the agreement of the numerical data (blue lines, made of 
data points) with the perturbative analysis (red lines) for small charge. Left: Merger between hairy 
black holes and RN AdS black holes along the line of marginal stability for 9 = 0.01. Right: Hairy 
black holes along a line of constant boundary condensate e = 0.1 for e = 3.2. 



infinity at fixed ( then we expect, on physical grounds that the limits 



9p{pX) 

fp{pX) 

Ap{pX) 

<Pp{pX) 



lim g{pL) 



fjpL) 
A{pL) 



lim 
lim 

L—^oo 



lim (t){pL) 

L— >oo 



(5.45) 



exist (see (3.8)) all appear to exist, so that the large charge hairy black hole may be rewritten 
in new coordinates r = Lp,t = T/L in the black brane form 



A 



-/p(p, C) dr^ + gpip, C) dp' + p' {dx' 
Mp,0 dr 



(5.46) 



We have numerical evidence that the limits (5.45) exist. We will leave the direct evidence for 
the next Section, where we will compare the planar limit of solitons and hairy black holes. For 
now, we verify the corresponding scaling of the mass with the charge and with the asymptotic 
scalar field. Working at the particular value e = 3.2, in Fig 21 (left) we present a graph of 

4 

M versus , where M represents the mass of the lightest black hole we have been able to 
construct at charge Q. In Fig. 21 (right) we present the plot of M verus the scalar vev, e, for 
the same solutions. These curves are both rather close to straight lines, demonstrating that 
the scaling described above works, at least for extremal black holes. 



Assuming that the limits (5.45) do exist in general, (5.46) yield the one parameter set of 



hairy black brane solutions labeled by C,. The entropy S of these solutions must scale like 
S = sL^, so that the ratio S/L^ has a good limit. These solutions should reduce to the black 
branes studied by HartnoU et al. [9,10 . We argued before that the same should be valid for the 
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Figure 20: Temperature vs. R (left) and Kretschmann invariant at the horizon vs. R (right) for 
e = 3.2 and a fixed e = 0.5. As the size of the black hole becomes smaller, the temperature decreases 
and it appears to become zero for R ~ 0.162. Simultaneously, the Kretschmann invariant seems to 
blow up as the black hole becomes smaller. Extrapolation of our data indicates that the blow up 
should occur for R « 0.154, which is close enough given the uncertainty. These results suggest that 
for 3 < < 32/3 hairy black holes posses an extremal singular limit. 



scaling limit of the second branch of the soliton; see Fig. 13 Indeed, the scaling of the mass, 
charge and scalar vev match roughly (notice that the large charges required to find the scaling 
become problematic for the numerics). 

6 Hairy black holes for > 32/3 



In this section we will study the solitonic and hairy black hole solutions of (1.1) in the parameter 
range > y. Let us first describe solitonic solutions. These solutions were constructed 
analytically at small charge (in a power series expansion in the charge) in [2]. In this paper 
we use numerical techniques to construct the solitonic solutions at all values of the charge. It 
appears to be the case that at all values of > y, solitonic solutions exist at arbitrary charge 
(i.e. they do not hit a Chandrashekhar instability at finite charge). Intuitively, for this range 
of parameters. Maxwell repulsion balances gravitational attraction without leading to collapse 
at the core. 

Let us now turn to the study of hairy black holes. In this range of parameters hairy black 
holes of charge q exist at all values of q. In |4||, the small q limit of these black holes was studied 
at fixed values of e^. In this section we supplement the study of j4] with a related but distinct 
perturbative construction of hairy black holes. We construct hairy black holes with charge ~ 6, 
in the theory with = y (1 + ^). In other words we construct hairy black holes whose charge 
is scaled to zero in a manner coordinated with the scaling of to y, in close analogy to our 
construction of hairy black holes in the previous section. 

It turns out that the nature of hairy black holes in this range of parameters undergoes a 
qualitative change when their charge exceeds a critical value. Below this critical value of charge 
qc2, the hairy black hole at its lowest mass ends on the soliton, as was predicted by |j4j. For 
charges greater than qc2, hairy black holes at lowest mass are extremal (similar to the behaviour 
of hairy black holes when 3 < < y). When = y (1 + ^) this critical value turns out to be 



39 




Figure 21: These plots M vs. Q (left) and M vs. e (right), for e = 3.2, represent our best estimate 
for where the extremal singular curve lies in the phase diagram, for large charges. In order to estimate 
the location of this curve, we take the points of smallest temperature for each line of constant e. The 
blue dots correspond to our numerical data and the red curves correspond to the best fit of our data 
to functions of the form M = A Q^^^ and M = B e'^ , respectively. We find that A ~ 1.0, B ~ 2.0 and 
C ~ 1.0. These results roughly agree with the scaling found in Fig. 13 for the second branch of the 
soliton. 



of order 6, and so is reliably captured by our new small 6 perturbation theory, though it lies 
outside the range of the perturbative expansion of |4]. 

The results of our new perturbative expansion, together with the results of |4] and certain 
other considerations, suggest a qualitative picture of the spectrum of hairy black holes at all 
values of q. In the rest of this section we then verify and fill in the details of this picture with 
the aid of extensive numerical simulations. 



6.1 Non-interacting model for small hairy black holes with = 

f(i+f) 

In this section, we predict the leading order thermodynamics of small hairy black holes for 

e' = f{l + e) (6.1) 

where ^ is a small positive parameter. 

We follow the methods of [5.1 to analyze these hairy black holes and define 

= aO 

fi' = l{l + aa9) (6.2) 
Qsoi = b9 

where a, a and b are positive numbers of order unity. In terms of these quantities, the thermo- 
dynamical expressions for the black hole and soliton are given by (modified versions of (5.4) 
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and (5.6)) 



fiBH = 2g J 1 + 



f^BH = 06* (2 + (1 + a)a9) , 
Qbh = 2aqc 1 + — — 



e + o{e^ 



(6.3) 



T 

Sbh 



(2 -a) 



271 



a9 , 



\3/2 



and 



rrisoi 





3b\ 


-qc + 




b 


262 




7 



(6.4) 



Note again that a < 2 (black holes with a > 2 formally have negative temperature and are 
unphysical; black holes with a = 2 have zero temperature and are extremal). 
As before, we consider systems with net charge 



and net masi 



T7 



m 



1 



q = q9, 



6 + {-q'^ + Am]e'^. 



Eq. (|5.9|) along with charge and mass conservations give the equations 

i]+o{e), 



a 
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Am 



a \7q, 
q = {2aq, + 6) + 0{e) , 
6 ah 862 
2gc qc 21 



(6.5) 
(6.6) 
(6.7) 



These equations can be solved for a and 6 

7 




{q + qc 



+ ^Am I + o{e) 



(6.8) 



Recall that the black hole component of the mix must have a < 2. From (6.5), this implies 
thatH 

h<^^{q + qc) + 0{e). (6.9) 



^^In other words, q and Am parametrize our hairy black holes. 

^*The inequality (6.91 is saturated when the hairy black hole is extremal. 
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It follows the negative root of (6.8) is the only acceptable solution. Positivity of a gives 

b<q + 0{e). 



(6.10) 



Inequalities (6.9) and (6.10) give rise to an important aspect of the phase diagram of these 
hairy black holes which we will now discuss. Two cases are possible: 



1. g < joiq + q, 



< 



7gc 



3 

In this case, the acceptable inequality for b is given by 

b<q. 

In particular, the above inequality is saturated at a = 0, i.e. the hairy black hole has a 
smooth soliton (a = 0) limit. This was predicted in In this case hairy black holes 
exist in the mass range (from (6.8)) 



o{e)>Am>-^-^-^ + oie). 

2qc 21 

The upper bound on Am is implied by the positivity of b. 



(6.11) 



2- g > Mq + qc 



The acceptable inequality for b is given by 



b<^{q + qc). 



This implies 



a > 



10 



7qc 



> 0. 



In particular the hairy black hole does not have a smooth solitonic (a = 0) limit. Note that 
the inequalities above are saturated when the hairy black hole is extremal (see footnote 



18 of the previous page). Hairy black holes exist when 



15 

At the lower bound the hairy black hole is extremal. 



(6.12) 



The complete phase diagram for this system at leading order is plotted in Fig. [22 
In summary, the non-interacting model for 
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^ yl + 6) predicts the following: 



1. Near extremal hairy black holes exist for all charges g > 0. 

2. For hairy black hole systems parametrized by q and Am (defined above) the soliton and 
black hole portions of the non-interacting mix have 



7 
10 



(q + q 



J_nq_7q, 
2qc \ 10 




(6.13) 



Am + 0(6) . 
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Figure 22: Range of allowed values of Am for hairy black holes with = y (1 + at leading 
order in small 9 (blue shaded). Hairy black holes exist for all charges g > 0. The upper limit, 
Am = denotes the onset of superradiant instabilities (red line). This is also the extremality 
line for pure RN AdS black holes. The green line denotes the soliton, Am = — ^ and 
the blue line denotes the extremal black hole. Am = —j^{q + 5'c)^- Note how the hairy black 
hole (at lowest mass) ends in the soliton for g < |gc and in the extremal hairy black hole for 
q > Iqc. 



3. Hairy black holes exist for mass ranges 

Oi9)>Am>-^-^ + 0{e), ifg<^; 

2qc 21 3 (g^^4^) 

0{e)>Am>-^iq + qcf + Oie), if g > ^ . 

15 3 

The lower limit of Am(g) undergoes a 'phase transition' at the critical value of charge 
qc2 = Below this critical value of charge, the hairy black hole at its lowest mass 
reduces to the soliton. Above this critical value of charge, the black hole component of 
the hairy black hole is extremal at its lowest allowed value of Am at any q. At this value 
of charge the lowest value of mass, Amc = Am(gc2) = "if- 
G.l.l Isotherms 



As we have explained in the previous section, nature of the lowest mass hairy black hole as 
a function of charge q undergoes a 'phase transition' at = |?c^- For smaller values of q 

the chemical potential of the solitonic solution is smaller than In this case the lowest 

mass hairy black hole is a soliton. Upon slightly raising the mass of these solutions at fixed 
q we nucleate a vanishingly small non extremal RN AdS black hole at the core of the soliton. 
The temperature of this black hole diverges in the zero size limit, so the solitonic solution is 
assigned temperature T = oo in the hairy black hole phase diagram. 



\/ i ^'^ chemical potential of the zero size extremal black hole 
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At g = the chemical potential of the soliton reaches the critical value y |. Upon raising 
the mass of the solution at fixed charge we nucleate a small extremal black hole at the center 
of the soliton. 

When q > qc^ the chemical potential of the soliton exceeds the critical value No zero 

size RN AdS black hole has a chemical potential larger than As a consequence, at these 
values of g, the soliton cannot be obtained as a limit of hairy black hole solutions. The lowest 
mass hairy black hole solutions in this limit are hairy extremal black holes, closely analogous 
to the solutions described in the previous section. These solutions are extremal and so have 
zero temperature. 

Note that the temperature of the lowest mass solution at charge q jumps from T = oo 
to T = as g increases past gc2- It is thus interesting to plot lines of constant temperature 
(isotherms) in the phase diagram for hairy black holes, in the neighbourhood of the lowest mass 
black hole at g = gc2- As T = and T = oo isotherms intersect at this point, it is natural to 
expect isotherms of every temperature to intersect at the same point. This is indeed the case 
as we now demonstrate. 

We define the parameter t as 

^ = a(2- a? = - {2a- — + 1 
Ve J ' a\ 7gc 

and look for curves of constant t. Using ( |6.13[ ), we can write 

20 (l20Am + 7 (4g (2g + ^6) + 3) ) 



7 (6y/6q + V2T^120Am + 7 (4g (2g + ^6) + 3) - 21 



+ 0(9). (6.15) 



Let us analyze the nature of this function in the neighbourhood of the critical point (|gc, ~f|)- 
Note that at this point the numerator and denominator both go to zero and hence t is inde- 
terminate. Slightly away from this point however, at g = |gc(l + ei). Am = — f|(l — ^2), t 
becomes (keeping upto linear terms in ei and €2) 

400(7ei + 5e2) 400c 



21(3ei + ^/20^/7^^~T5^) 21(3ei + V20c} 
where c = 7ei + 5e2. If we choose ei = — |-\/20c + 7c and take c — )■ 0, t reduces to 



(6.16) 



The analysis above clearly demonstrates that it is possible to approach the critical point 
(|gc, — f|) via some path so as to attain any temperature t as is required (by an appropri- 
ate choice of 7). In other words, all isotherms of this system must pass through this point. 



This is more clearly seen in Fig. 23 



6.2 Perturbative analysis of hairy black holes 

We have followed the procedure described in Appendix |B. 4 to construct the hairy black hole for 
= y (1 -|- ^) in a perturbative expansion in 9. We have listed the explicit results of our con- 
struction in Appendix [Dl With these solutions in hand the evaluation of their thermodynamic 
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Am 




Figure 23: Isotherms for hairy black holes at leading order. These have been plotted for 
various values of t as indicated on the graph. The black curve is the extremal black hole 
Am = — -^{q + Qc]"^ ■ The green curve is the soliton line Am = — Hairy black holes 
exist in the pink shaded region. The two insets show zoomed in regions of the plot. Note that 
for q < l^qc, the isotherms lie above the soliton, and approach it as t — )■ oo. For q > l^qc, the 
isotherms line between the soliton {t = oo limit) and extremal hairy black hole [t = limit). 
In particular, all the isotherms intersect the soliton and extremal hairy black hole curve at 
g= |g„ Am = -f|. 
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charges and potentials is straightforward. Using formulas given in ^5.2, we find 



M 

Q 

T 



e 2a 



8k 



TT /3 
4V2 



9 

e{2a + 



1024afc 



a + 



9232fc^ 



189 59535 
8k\ 9^ (322560aA; - 59535a + 4192P + 26460fc) 



9 



59535 



+ O (6^) , 



3 {8k -21)9 e2(^11088(1642a + 525)A; + 4584195(2a + 1)2 + 135136P) ^.^3. 
2 14v^ 4074840^6 ^ ^ 

y^(42a - 8A; + 21) 
427rv^ 

^3/2 (-5544(2291a + 945)A; - 1528065(3a(a + 1) + 1) + 965368^) 



+ 



30561307rv/^ 



(6.18) 



We have verified that these quantities satisfy the first law of thermodynamics 

dM = TdS + jjdQ . 



To verify the results of §6.1[ we work with rescaled and shifted mass Am and rescaled charge 



9a + Ak 



Am 



3^6 
4 



(6.19) 
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A;(126a+ 16A; + 63) + C(^). 



As before, identifying b = 7^ (see ^5.2.4), the formulae (6.19) reduce to 



Am 



q = {2aq, + h) + 0{9), 
h ah 862 
2gc Qc 21 



+ 0{9). 



(6.20) 



These are exactly the formulae (6.6) and (6.7). The hairy black hole can therefore be considered 
as a non-interacting mix of an RN AdS black hole and a soliton. 



Breakdown of perturbation theory at extremality For q < qc^ our perturbative con- 
struction of hairy black holes applies at every value of the charge. When q> qc2i however, the 
value of the scalar field blows up at the horizon of the extremal hairy black hole, which indicates 
a breakdown in our perturbative expansion near extremality, exactly as in the previous section. 
Also as in the previous section our perturbative construction (and so perturbative predictions 
for thermodynamics) is reliable down to exponentially low temperatures T ~ q-^I^"-^ , 



6.3 Hairy black holes at = 



3 



As we have explained above, when e = ^{1 — 9) with 9 small and positive, (|l.l|) hosts hairy 



black holes of charge greater than or equal to qc9. The hairy black hole of the lowest mass at 
any given charge has vanishing temperature. It may roughly be thought of as a small extremal 
RN AdS black hole embedded in a soliton. 
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On the other hand when = y(l + 9) with 9 small and positive, (1.1) hosts hairy black 
holes at all values of the charge. The nature of lowest mass hairy black holes, at any given value 
of the charge, however changes discontinuously at charge equal to T^qc9. Below this charge, the 
lowest mass hairy black hole is a regular soliton and has infinite temperature. Above this charge 
the lowest mass hairy black hole is at zero temperature and may roughly be thought of as a 
small extremal RN AdS black hole embedded inside the soliton. 

What is the situation at = y? This special value of the charge may be reached as the 
6' — )■ limit of either = y (1 — 6') or = y (1 + either limit suggests that we should find 
hairy black holes at all values of the charge, and that the lowest mass hairy black hole at any 
given charge should be a zero temperature solution that may roughly be thought of as a small 
extremal RN AdS black hole embedded inside the soliton. In Appendix [C] we have verified that 
this picture is correct, by explicitly constructing all hairy black holes at small charge in the 
theory with = y. 

6.4 ^C2(e^) larger values of e 

We have used our perturbative expansion above to demonstrate that 

7 /3e2 \ 9127gc /3e 




We^)^-,.^--lJ-^^--lJ H-OI 1^-11 |. (6.21) 

It is, moreover, an interesting fact that Q'cale^) is related to a simple property of the solitonic 



solutions of (1.1), and so may be computed quite easily even at larger values of using the 
following considerations. When q < Q'cale^) the solitonic solution is a limit of hairy black hole 
solutions. In other words, at these charges there exist hairy black holes that may be thought 
of as an infinitesimal RN AdS black hole sitting at the center of the soliton. As the RN AdS 
black hole is infinitesimal, the resultant solution may accurately be thought of as a mix of the 
soliton and the RN AdS black hole at every value of the charge q (not just at small q). As 
we have explained in earlier sections, the conditions for stationarity of such a mix are that 
the chemical potential of the small RN AdS black hole matches the chemical potential of the 
soliton. However the chemical potential of infinitesimal RN AdS black holes is bounded from 

above by It follows that if g < then the chemical potential of the soliton at that charge 
is less than -i /|. The charge qc^ is determined by the requirement that the chemical potential 



of the soliton at charge precisely equals ^J~^. 



6.5 Numerical results 



The set-up for constructing hairy black holes numerically was already detailed in subsection 5.3 
Here, we present the results for the range > y . Recall that, for this range, the soliton familj 
of solutions exists for all charges. 

The most remarkable feature of the phase diagram, which we represent in Fig. |24| is the 
lower mass bound of the hairy black holes. The Fig. confirms the phase transition predicted 



from the perturbative analysis. It should be compared to Fig. 22 there, close to y was 



used, while we now have e = 4. For charges smaller than a critical value, the hairy black 
holes have a zero size limit which is the soliton branch. For higher charges, however, the lower 
mass bound consists of extremal but singular solutions, as happened already for 3 < < y. 



This behaviour, which should occur for arbitrarily large charges, was described by Ref. 29 
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Figure 24: Left: AM vs. Q for e = 4. The black curve corresponds to the sohton and the blue region 
corresponds to hairy black holes. The red curve is the line of marginal modes of the linear problem 
discussed in section 4^, and the shaded region is the area occupied by RN AdS black holes. The black 
curve now lies entirely below the region occupied by RN AdS black holes and it continues for arbitrarily 
large values of Q. The blue curves end at the black curve for small charges {Q < 0.75), indicating 
that the soliton branch is the endpoint of the hairy black holes family. However, for large charges 
{Q ^ 0.75), the lower mass bound of the hairy black holes is below the soliton. These numerical results 



reproduce the behaviour predicted in Fig. 22, now for e = 4. Right: same data with the entropy as 
third axis. The red surface is the strip of RN AdS black holes between the marginal line of stability 
and extremality. 



in the planar case. It was found there too that the extremal solution is singular. In order 
to make a quantitative comparison with the perturbative prediction of (6.18), we represent in 
Fig. 25 hairy black holes along the merger line (upper mass bound) and along a line of constant 
boundary condensate e = 0.1, for values of close to y. The agreement for small charge is 
clear. 

We perform another check of the numerical results in Fig. 
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choosing e 



3.4, close to y. 

On the left, we verify that the soliton and the merger (between RN AdS and hairy black holes) 
curves agree for small charge with the perturbative predictions of (6.14) and (6.18). We also 
want to check that curves of hairy black holes whose limit is a singular extremal black hole 
rather than the soliton, i.e. terminating for Q > Qc^ predicted in (6.21), cross the soliton line 
and have a zero-temperature endpoint; see right plot. This seems consistent for e = 0.7, and 
the opposite situation (endpoint is the soliton and temperature diverges) seems consistent for 
e = 0.2 and e = 0.3. Unfortunately, for the curves closer to the transition, e = 0.4 and e = 0.5, 
the results are inconclusive. The intricate behaviour near the transition presented a technical 
obstacle in our effort to approach numerically the endpoints of the curves. 

As in the previous section, we expect that hairy black holes at > y admit the planar 



scaling limit (5.45). We provide numerical evidence that such a limit exists for e = 4 in Fig. 27 

4 

In particular, we confirm that M scales as and, equivalently, as e for large charges. These 



results agree with the soliton scaling, according to Fig. 10 



We can do better and contruct approximately the one parameter set of planar solutions 
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Figure 25: These plots M vs. Q represent the agreement of the numerical data (blue lines, made of 
data points) with the perturbative analysis (red lines) for small charge. Left: Merger between hairy 
black holes and RN AdS black holes along the line of marginal stability for 9 = 0.01. Right: Hairy 
black holes along a line of constant boundary condensate e = 0.1 for e = 3.33. 



(5.46) through an analysis for large charges. Our results are presented in Fig. 28 for e = 4. We 



show the behaviour of the metric function fp for fixed ( (here ( 
e = 4, according to the scaling in Fig 



27). 



0.8, which corresponds to 
We present the results for both the solitons and 
the extremal hairy black holes, which we approximate by the lowest mass black holes we can 
find numerically for a given charge. The two limits appear to coincide, within our accuracy. 
This provides evidence that the solitons and the extremal hairy black holes have the same large 
charge limit, which we expect to be the zero temperature limit of the black branes studied by 
HartnoU et al. [9,10 . We expect also that our non-extremal hairy black holes have as a planar 
limit the black branes at finite temperature. 
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Figure 26: Left: AM vs. Q for e = 3.4. The black curve corresponds to the sohton, which agrees 
for small charge with the perturbative prediction represented by the green dashed line. The red curve 



corresponds to the line of marginal modes of the linear problem discussed in section 4.2, which agrees 



for small charge with the perturbative prediction represented by the magenta dashed line. The other 
lines correspond to hairy black holes at different fixed values of the boundary condensate e (see legend 
on the right plot). represents the perturbative prediction for the change in the lower mass bound 
of hairy black holes. It is clear that, for e = 0.7, the hairy black holes extend to the the left of the 
soliton curve; the same should happen for e = 0.4 and e = 0.5, but it was not possible to reach a 
sufficiently small value of the charge with our numerical approach. Right: T vs. R for e = 3.4. The 
curves correspond to hairy black holes at the same fixed values of e. The temperature should diverge 
for e = 0.2 and e = 0.3 as R vanishes, and it should vanish in the other cases at a small but non-zero 
value of R. The numerical results are inconclusive because there is an intricate behaviour around the 
transition, so that we cannot reach sufficiently small values of R with our technique. Nevertheless, 
that does seem to be the case at least for e = 0.5 and e = 0.7; for e = 0.3 the curve might oscillate 
but end up turning down at smaller R. 



A Non-existence of regular planar solitons 

In this appendix we will provide an argument for the non-existence of regular solitons in the 
planar limit. 

In the planar case, the ansatz for the metric, gauge field and scalar field is 

ds^ = -f(r)dt^ + — + r^dx'?, Aa = A(r)(dt)a = 0(r) G M (A.l) 
9{r) 

Following the same procedure as in the main text we can derive 2nd order (quasi-linear) 
equations of motion for 0(r), A{r) and 0(r). The general structure of the equations is as follows: 

Q" + M^^^^'\Q)rQ',{Q'f + M^,^^'\Q)Q',Q[ + M^^^/{Q)^Q'^ + _Lm(4).(Q) = 0, (A.2) 

where the objects M(q, £ = 1, . . . , 4, are non-hnear functions of the variables Qi{= f, A, 0) and 
M(i) and M(2) need not have any special symmetry in the pair of indices {kl). 
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Figure 27: These plots M vs. Q (left) and M vs. e (right) for e = 4, represent our best estimate for 
where the extremal singular curve lies in the phase diagram, for large charges. In order to estimate 
the shape of this singular curve in the phase diagram, we take the points of smallest temperature for 
each line of constant e. The blue dots correspond to our numerical data and the red curves correspond 
to the best fit of our data to functions of the form M = AQ^I'^ and M = B ^ respectively. We find 
that A ~ 0.8, B C ^ 1.0. These results agree with the scaling found in Fig. 10 for the soliton. 



The idea is to solve the equations near the origin r = 0, assuming that the functions Qi 
have an expansion in non negative integer powers of r: 



20 



N 



.N+l\ 



(A.3) 



1=0 



for r — )■ and where the qi^e are constants that are determined from the expansion, i.e., the 
values of the functions and their first derivatives at r = 0. By an induction argument then we 
will show that in fact the expansion above is in even powers of r only. 

Indeed, solving (A. 2) near r = we find that the first few terms in the expansions contain 
only even powers of r: 



(r) = - + 0(r2), A(r) = 0(r2), f(r)=0{r^). 



(A.4) 



As discussed in 29 , one can use the scaling symmetries of theory and the various fields in the 
planar case to set the leading terms in A and / to 1. Therefore, we have: 

2 e 



r) = r^ + 0(r'), A(r) = + 0(r^) , f(r) = r^ + 0(r^) 

e 6 



(A.5) 



Now we assume that up to ©(r^"') the various functions have an expansion in even powers of 
r, but at the next order there is an odd power of r: 



2n+2\ 



(A.6) 



^"Non-integer powers of r would lead to singularities in the derivatives of the metric, and therefore in the 
curvature tensor (or derivatives of it). 
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Figure 28: Metric function fp of a planar solution (5.46), constructed numerically from the scaling 
(5.45), for e = 4, i.e. C, = 0.8. The black line represents the scaling solution of the solitons, and the 



green dots represent the scaling solution of the extremal hairy black holes. 

Next we have to show that the ©(r^""*"^) term above vanishes because of the equations of motion 



(A.2). 



From (A.2) and the fact that, by assumption, we have determined the solution up to 0{r 



2n\ 



we see that the leading term appears at 0(r^" ^). Clearly the terms in (A.2) with M(i) and 



M(2) cannot contribute at this order: either they give terms which are of higher order or terms 
which are even in r. Therefore, the only contribution at this order can come from Q'l and the 
last two terms. Note that the combination ^ Q'^ gives terms with positive even powers of r up 

to 0(r^"^^). Depending on the form of Mf^^-^r' we may also get a contribution from these other 
terms. Using the explicit expressions, one finds that the 0(r^"~^) contribution of these terms 
is 

i = (j): 3(2n + 1) q^^-zn+i + ^[^ Um+i + - 3)g/,2n+i] , 
i = A: [3{2n + 1) + 4] gA,2n+i + 2(2n - l)g/,2n+i , 
i = f : (2n + l)g/,2n+i • 
To close the argument we need the explicit form of M(4). The components are: 

3e^A'^(j) 6e^Af(j)'^ 12/ (2/ + e^V) 



(A.7) 



(4)0 



Ml 



M, 



6/ + 6^202 ' ^ 6/ + 62^202 ' (4) / ^ ^2^202 

One can see that the 0{r'^^~^) contribution that these terms give rise to is: 



(A.8) 



M, 



(2n.-l) 



(2n-l) 



M, 



(2n-l) 



-4g/,2n+i. (A.9) 



'(4)0 — e(2gA,2n+l '?/,2n+l) , ^"(4)A ~ -^^A,-Zn+i, ^"(4)/ 

Therefore, at 0{r'^"'~^) we are led to a homogeneous system of linear equations for the 

Q'i,2n+l's: 

{2n + l)2n g0,2n+i + 3(2n + 1) g0,2n+i + 2 QA,2n+i + (2?^ ^ 3)g/,2n+i] + e(2 qA,2n+i - g/,2n+i) = , 

(2n + l)2n qA,2n+i + 3(2r2 + 1) g0,2„+i + [3{2n + 1) + 4] gA,2n+i + 2(2n - l)g/,2n+i - 4 qA,2n+i = , 
(2n + l)2n g/,2n+i + {2n + l)g/,2n+i - 4 g/,2n+i = . 

(A. 10) 
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But clearly this system is non-degenerate and the only solution is the trivial one: g<^,2n+i = 
(lA,2n+i = Qf,2n+i = 0. Therefore, by induction, we have shown that the near r = expansion 
can only contain even powers of r if we demand regularity. 

Now we can integrate out the equations of motion and we find that all coefficients in this 
expansion are be determined in terms of e only (we have constructed such an expansion up to 
0(r^^). Therefore we see that this near horizon expansion has no free parameters. Ultimately 
this is a consequence of the regularity requirements that we have imposed. If there was a solution 
in the whole spacetime, then this near "horizon" solution should be matched to a far region 
solution with the correct asymptotic behavior. However, one can easily check (numerically) 
that this is only possible at isolated values of e. This is precisely the result that [29] found: 
only at one particular value of the scalar charge there is regular soliton. Therefore, we conclude 
that for generic values of the scalar charge e there is no regular planar soliton. 



B Perturbative construction of small hairy black holes 

for e2 = f (1 - 9) 

In this section, we present a detailed description of the perturbative procedure used to construct 
small charge hairy black holes in an expansion in 6. We also present a detailed listing of our 
findings. 

B.l Setting up the perturbation theory 



As discussed in §5.2.3 RN AdS black holes form the starting point of our perturbative expan- 



sion. In the following paragraph, we briefiy review RN AdS black holes. 



RN AdS Black Holes The AdS-Reissner-Nordstrom black holes constitute a very well 
known two parameter set of solutions to the equations of ^5.2.1 These solutions are given 
by 



V{r) 



A{r) 
0(r) 



-V{r)dt'^ 



+ r^d^ll 



1 + R^ + -/i' 



2 
3 



1 - 



21 r 



1 + r^ + R^ - 



2/i2/?2 



(B.l) 







R^ 



where is the chemical potential of the RN AdS black hole. The function V{r) in (B.l ) vanishes 
at r = R and consequently this solution has a horizon at r = R. We will require that R is the 
outer horizon of our solution; this imposes the restriction (see El) 



^'<1{1 + 2R^) 



(B.2) 
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In this paper we are interested in small, near extremal RN AdS black holes. Following ^ 5.1.4 

we set 
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The RN AdS black hole becomes 

ds^-- 



= a9 
_ 3 
" 2 ' 



1 + aa6 



(B.3) 



V{r) 
A{r) 



V{r)de + — - + rHnl 
^ ' V{r) 

a9{aa9 + 1) 



+ a9 + r'^ 



(B.4) 



[1 + aai 



aO 



The condition (B.2) reduces to 



a<2 



(B.5) 



The thermodynamics of these RN AdS black holes was reviewed in ^5.1.4 



As we have explained in the main text (see ^5.2.3), our perturbative construction of small 
hairy black holes uses (B.4) as the starting point. Our construction proceeds by systematically 
correcting the background solution (B.4) in a power series expansion in \/9 in three different 
regions (near field, intermediate field and far field) and finally patching the full solution together 
by matching in regions of overlap. We now describe our procedure in more detail in each of the 
three regions. 

B.1.1 The far field region r \fQ 



In this region the starting RN AdS solution (B.4) is a small perturbation about global AdSs 
space. The scalar condensate - of order a/^ - is also a small perturbation about global AdSs 
space. Consequently we seek solutions of the form 

oo 

n=0 

oo 



n=0 

oo 



out I 



n „out I 



(B.6) 



n=0 



out 

'(27i+l)/2 



n=0 



As the starting point of our perturbation theory is the RN AdS black hole, we have 



aOiaaO + 11 



+ a9 + r^ + 1 



Mr) 



(B.7) 



'1 + aal 



1 - 



a9 
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The superscript out emphasises that this expansion is good at large r. Notice that /g''*, 
and Aq"* depend on two parameters a and a. While a is one of the parameters of the hairy 
black hole solution that we wish to find, a is not. Our perturbative procedure will determine 



a as a function of a and k (see (5.32)) in a power series expansion in 9 
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a = an{a, k)9"' 



fB.81 



ra=0 



We pause to explain a slightly confusing aspect of the notation of (|B.7|). As is apparent from 
(B.7), the starting point of our perturbative expansion 



fout 

Jo ' 



^out 



and A'^*' are each given 



functions of 6', each of which starts out at Oil) and admits a regular expansion in small 6. The 
higher order corrections, however, /°"*, (7°"* and ^4°"* {n > 0) in the expansion (B.6) are each 



taken to be independent of 9. We adopt this convention to east matching with the intermediate 
field region, 
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The far field perturbation expansion is expected to break down when the base (zero order) 
solution deviates substantially from the metric of AdSs space. A cursory inspection of ( |B.7 ) 
reveals that this happens at r ~ \/a9. 



B.1.2 The intermediate field region r ^ 1 and r — R — Rin ~ 



^3 



The intermediate region is defined by C,y a9 <^ r — \'a9 <^ 1 (see (B.3)) . Over these length 
scales the small black hole is far from a small perturbation about AdSs space. Instead the 
simplification in this region stems from the fact that we focus on distances of order 0(^/9). 
Over these small length scales the background gauge field, which is of order unity, is negligible 
compared to the derivatives which are of order 4=. A second simplification results from the 
fact that we insist that r — R^ R — Rin, i.e. we do not let our length scales become too small. 
At these distances the black hole that we perturb around are effectively extremal at leading 
order. Moreover the black hole may also be thought of (at leading order) as a small black hole 
in flat rather than global AdS space. 

In this region it is convenient to work in a rescaled radial coordinate y = -f^ and a rescaled 



time coordinate r = -y=^. Note that the intermediate field region consists of space time points 



with y of order 0{1). Points with y of order (or larger) and y — 1 of 0{9) (or smaller) are 
excluded from the considerations of this subsection. 



^"'^Our perturbative expansion has an ambiguity, corresponding to a ^ dependent redefinition of a in the RN 
AdS solution used as the starting point of the perturbative expansion. This ambiguity of course has no physical 
significance: it corresponds to different partitionings of the same final solution into the 'zero order' part and the 
'perturbation'. In this paper we fix this ambiguity by requiring that the correction to the gauge field at infinity 
vanishes at every nonzero order in perturbation theory. This convention removes the ambiguity in perturbation 
theory, and gives physical significance to a. Concretely, the chemical potential of our solution is given in terms 



of a via (B.3) . 

^^This convention differs from the usual one (that all /^"*, 5™* and are independent of 9) by a simply 
shift of all higher order terms in the 6 expansion. 
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The metric and the gauge field of the background black hole take the form 



ds^ = aO 



V{y) 
AM 



-V{y)dr^ + ^ + yHnl 



y^ 



y^ 



3 /I 

-a9{l + aa9) 1-^ 
2 V 1/ 



(B.9) 



a 



n=0 



Note the metric and gauge field in (B.9) are not independent of 9 even at leading order 



in small 9. However the leading order dependence of the metric and gauge field on 9 is very 
simple; g^^ ~ d{9) and ~ 0{-\f9). Now the action for our system 



S 



SttCs 

can be rewritten in terms of the rescaled variables 

I I mid 



(B.IO) 



(B.ll) 



as 



D 



S 



mid 



02 

SttG, 



(Tx 



^Qprnid 



-detig^i'^) 



If 

4 i'" 



mid 



mid\fj,u 



+ \D 



mid 1 mid 1 2 



+ Q9 



(B.12) 



The net effect of the rescaling is three fold. First we multiply the action of the system by a 
constant factor, which does not affect classical equations in any way. Second the cosmological 
constant term is multiplied by a factor of 9, implying that the cosmological constant is a small 
perturbation to the dynamics in the intermediate field region. Third the coupling of the Maxwell 
term in the covariant derivative is multiplied by a factor of v^, implying that the scalar and 
gauge field are almost decoupled in this region. 

Now f?^**^ and A^*"^ have a good ^ — )■ limit on the background RN AdS solution. This 



suggests that 5'™'^ and A^^'^ and should admit a standard expansion in the amplitude of the 
scalar field in the intermediate field region, when expressed as functions of the rescaled variable 
y. In other words if we define 



ds^ 

j^mid 



a9 [-r\y)dT' + g"''\y)dy' + y'dnj] 



(B.13) 
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then we can expand 



n=0 
oo 

r (B.14) 



n=0 

oo 

Lmidi 



where 



n=0 



/o(y) (B.15) 
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The perturbative procedure that determines all the higher order terms in (B.14) is now just 
a standard expansion in the amplitude of the scalar field 0, very much like in the far field 
region. After solving for the hairy black hole in the far and intermediate field regions we will 
need to match the two solutions. In order to facilitate this matching it is useful to rewrite the 
intermediate field metric and gauge field in terms of the far field coordinates. We find 



ds' = -r'\r/R)dt' + g"'"'{r/R)dr' + r'd^i 
At = A'^'^r/R), Ar = A = 0, 0(r) = 0™'^(r/i?) 



(B.16) 



Comparing with (3.1) we see that we must directly match the functions f°^^{r) with J™**^ (r/R), 
g°^\r) with c/™"^ (r/i?), A°"*(r) with A'^''^ {r / R) and (p°^\r) with 0™'^ (r//2). 

Finally let us examine when we expect the intermediate field expansion to break down. We 
expect a break down when y"^ is order ^ or greater, as the last term in the expression for Viy) 



in (B.9) cannot be treated as a small perturbation in this region. Moreover the two roots of the 
equation Viy) = are given hj y = 1 and y = 1 — ^a(9 for y — 1 ^ ^a(9 the two roots appear 
coincident, and dynamics effectively takes place in the background of an extremal black hole. 
At ?/ — 1 ~ ^cl(0, however, the equations see the deviation from extremality in an important 
manner, and the intermediate field expansion is expected to break down. 



B.1.3 The near field region r — R <^ R 



In this subsection we now turn to the near field region r — yaO <C va6 (see (B.3). We will 



work in terms of a further rescaled radial coordinate z = ^-4- Note that the black hole horizon 



^^As in the previous subsection /™*'^, g™''^ and A™'^ have a specified dependence on 9, whereas all other 
higher order terms in the perturbation theory /^"^, and A^*'' (n > 0) do not depend on 9. As before, this 
convention eases matching the intermediate field regions to the far-field and near-field regions. 
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occurs at 2; = and that points with z of order unity are located at r — i? ~ i? — Rin ~ 
C-R^ ~ or y — 1 ~ QW" ~ aQQ. It is also convenient to work with the new time coordinate 
T = \faBC^t = a(9T. As in the previous subsection, the background gauge field makes a small 
direct contribution to dynamics in this region. However deviation of the black hole metric 
from extremahty (and the difference between an AdS and flat space black hole metric) are all 
important in this region, and have to be dealt with exactly rather than perturbatively. Let us 
now see all of this in more detail. 

In the new coordinates, the RN AdS metric and gauge field take the form 



ds^ = a9 



-V{z)dT' + + (1 + aCez^dQl 



^ ^ 0,2(25)2 \^ {azCO + l)^J V V / {l + az(9)^J 

= 2z{l + 2z) + 0{e) 
00 



a — 



As in the previous subsubsection the metric and gauge field do not directly have a good ^ — > 
limit. However the rescaled gauge field and metric 

„ n^near 

= \feAl^''^ (B.18) 

(j) = 0"^"^ 

are well defined in the small 6 limit; moreover exactly as in the previous section, the rescaled 

variables obey the same equations of motion as in rescaled quantities apart from a rescaling of 
the cosmo logical constant (which sets it to 0{6)). It follows that the functions defined by the 
expansions 

ds"^ = ae [-P'"''-{z)dT'^ + g'""''-{z)dz'^ + (1 + aCezfdnf\ 

At = V^A'^^'^'iz), A, = A = 0, = (^""""^(^z) ^^'"^^^ 

should admit a power series expansion in \fd, 

00 

n=0 

00 

V (B-20) 



n=0 

00 

0--(z) =V^Er0^i:;i)/2(^) 



n=0 
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where 



a^CW \ {azCe + 1)V V (1 + 

1 



/r"^'(^) (B.21) 



CV2 a202 (l + aC^^^)' 

0r'(z) =0 



The perturbative procedure that recursively determines higher order terms in (B.20) is now 
simply a standard expansion in the amplitude of the scalar field (p, just as in the intermediate 
and far field regions. 

In order to match with the functions gmid^ j^mid ^^'^^ it is useful to rewrite the 



configuration (B.19) in the original coordinates r and t. We find 



We conclude that, for matching purposes 



(B.22) 



J out 


rmid 
~/ ~ 




gOut 




gUear 




A mid 




.out 


imid 


inear 



(B.23) 



In order to facilitate this matching we simply define further rescaled functions 

r{z) ={acef r'^'Xz) 

^ aX^O^ (B.24) 

so that (for matching purposes) 

yo«t ^ jmid ^ pn^ gOut ^ ginid ^ gin ^ j^out ^ j^mid ^ j^in ^out ^ ^mid ^ ^in 

Below we will find it convenient to use the functions in actually implementing our pertur- 



bative expansion. It follows from (B.20) and (B.24) that and A*" admit the expansions 
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24 



n=0 

"=° (B.25) 

n=0 



n=0 



where 



a2^2 \^ ^azCe + 1)V V (1 + a-^C^)^ 

(B.26) 



B.1.4 Subtleties in the near field region 

There are two related further subtleties in working out the perturbative expansion in the near 
region. First recall that the radial coordinate r employed in this paper has geometrical sig- 
nificance; it parametrizes the volume of the surrounding at that point. For this reason 
re-parametrizations of r do not, in general form a symmetry of the equations in this paper. At 
leading order in the near field region, however, the metric metric takes the form 



(B.27) 



Note in particular that the size of three sphere (at leading order) is a constant independent of 
z. Therefore, at this order, coordinate re-parametrizations of z are a symmetry of the problem. 

Now consider working out the equations of motion for the near field metric and gauge field 
perturbatively in 6. At leading order in this expansion, it follows from our discussion above 
that the equations in the near field region will admit a whole functions worth (instead of 4 
numbers worth) of solutions, parametrized by any 0{6) redefinition of z coordinate. In order 
for this to work out in practice, it must clearly be the case that the differential equations of 



^5.1.4 are not all independent at leading order in 9. Indeed the relation that we expect between 



the equations is easy to deduce. A coordinate redefinition of z 

z = z' + 6{z') 

induces the following variable changes 

5gzz = -Siz)d^gzz{z) + 25'{z)gzz{z) 

SgTT = -Siz)dzgTT (B.28) 
SAt = -6{z)dzAT 
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The appropriate factor of ( has been absorbed into the functions f^{z), g^{z) and A^j!^{z) 



60 



The change in the action induced by these variable changes must vanish (at leading order in 
6) because it is effectively a coordinate transformation. Independently, however, the change in 
the action due to these field variations is given by 

6S oc E^^dgTT + E''6g,, + M^SAt = 

where E^^, E'''' and are the TT, zz of the Einstein Equation and the T component of the 



Maxwell Equation. Using (B.28), this gives 



29. iy^-E,:) - — Mr = ^ [Err + (V^n E..) (B.29) 



where Vq'^{z) = 2z(l + 2z). The net upshot of this discussion is that we expect (B.29) to be 
an identity at leading order in 6', and it is not difficult to directly verify (see below) that that 
is indeed the case. 

The equations admit a second subtlety of similar nature. It turns out that to leading order 
in 9, the action is a function only of the square root of the determinant of the metric in the 



z and T directions, rather than separately a function of qtt and Qzz "^ For this reason metric 
variations of the form 

leave the action unchanged (at leading order). Consequently it must be that 

Ett + {Vn'Ezz = (B.30) 

is an identity at leading order, and this is easy to directly verify. 

Thus we have deduced, without even doing any calculations that two linear combinations of 
the three functions in the metric and the gauge field should be undetermined at leading order 
in the perturbative expansion. 

How do we deal with these subtleties. It is very important that the subtleties referred to 
above are purely leading order phenomena; at higher orders in 6 coordinate redefinitions in z 
are not symmetries of our system because the size of the sphere is not constant at subleading 
orders in 6; furthermore at higher orders in 6 the action depends separately on qtt and Qzz 
instead of only on their product. 

In order to see the import of all of this, consider perturbation theory at (9(6'^). As the 
homogeneous part of the equations are same at every order, the same linear combinations 



(B.29) and (B.30) of the (9(^^) fluctuations disappears from (i.e. is undetermined by) the 
second order equations. However, the two functions that were left undetermined in the 0{6) 
calculations now appear in the 0{9^) calculations. 

More precisely, at every order other than the leading, we actually do have as many equations 
as variables. The variables, however, consist of one unknown function at that order coupled 



with two unknown functions of the previous order It turns out that the perturbative near 



^^In order to sec this, we note that the metric at leading order is simply AdS2 x S'3. Therefore, the term 
TZ[g] in the metric simply topological. Further, we will later see that the 0{Vd) calculations at near field set 
^1/2 (■^) ~ ^- Consequently, the only non-trivial part of the Lagrangian at leading order 0{9) is ^/gF^^F'^'^ = 
\/99^'^ 9^'^' {FtzY , which is a function only of ^ 

^^These are the unknown coordinate transformations and unknown ratio gxr/dzz 



61 



field equations at order 0{6"^) take the form 



2^ 

3dz^ 



9m{z) + 



V^^iz)^Am^,iz) 



wliere 



1 d 

3CV]^d^ 

2('^z{l + 2z) 



^m-l(^) 
^m-l(^) 

TKz) 



(B.31) 



Terms on tlie RHS of the equation above are referred to as source terms. The source terms 
T^_i{z), tI^_^{z) are determined completely by An-i, fn-i and gn for n < m. On the other 
hand he source term T^{z) depends on An, fn and Qn for n < m. In order to implement 
perturbation theory, we solve the first two equations in (B.31) at m = 2 along with the last 



equation in (B.31) at m = 1 (this determines Ai, fi and gi). At next order we solve the first 
two equations in (B.31) at m = 3 along with the last equation in (B.31) at m = 2. This 



procedure continues recursively. 

Finally, let us examine when we expect the near field expansion to break down. We expect 
a break down when z is order | or greater, i.e. r — y/aO ~ y/aO, as the last two terms in V{z) 



in (B.17) cannot be treated a small perturbation in this region. 



B.1.5 Summary 

According to the discussion above, in order to evaluate the hairy black hole solution we must 
implement the perturbative procedure discussed above separately in each of the three regions 



(near, intermediate and far), implement boundary conditions of normalisability at infinity (5.31 ) 
in the far field region and regularity at horizon in the near field region, and then match all the 
solutions into a smooth whole. We will now implement this procedure in detail upto order 6. 

B.2 Perturbation theory at 0{Ve) 

In this subsection we work to order 0{V6). The procedure described in this subsection applies 

with minor modifications to the perturbative construction at 0{\fQ ^ 
equation at 0{V6) is 



). The only non-trivial 







where = V^^ — ieA^j^ is the linearised gauge covariantised Laplace equation about the back- 
ground (B.l). We will now solve this equation subject to the constraints of normalisability of 



infinity, regularity at the horizon, and the requirement that 



+ 0{l/r^) at large r. 



B.2.1 Far field region r ^ 

Let us first focus on the region r ^ \/6. In this region the background (B.4) is a small 



perturbation about global AdS space with At = \J\. A stationary linearised fluctuation about 
this background is gauge equivalent to a linearised fluctuation with time dependence 



^^The time dependence is e 
leading order in 



-iJzl2JzilZt 



-4it 



where we have used = ^(1 — 6') and worked to 
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about global AdS space with At = 0. The required solution is simply the ground state excitation 
of a massless minimally coupled scalar field about global AdS. The far- field scalar equation at 

o{Ve) is 



d 



d 



dr V (1 + ^^)^ dr 
The normalizable solution is 



Vk 







1 + r 



2\2 



(B.32) 
(B.33) 



the overall normalisation of the mode is set by the requirement (5.32) 



+ O (l/r^) 



This completes our perturbation analysis at 0{V6). 

The far-field scalar equation at O ^ ^ generalises to 



d 



d 



dr V (1 + '^^)^ dr 



[(1 + 



,2\2iout 



^{2m+l)/2\ 



-P(2m+l)/2i^' 



(B.34) 



The source term P; 



out 



(2m+l)/2 



V) is determined by the lower order solutions. We could therefore 



iterate this process to generate 0°2m+i)/2(^) ^^^^ desired m. 



It turns out that the expressions 0('2m+i)/2(^) increasingly singular as r — )• and the 
most singular piece scales like ^2{^-i) {rn > 0). In other words, the expansion of 0°"* is really 
an expansion in ^ (upto log corrections and an overall factor). The expansion therefore breaks 
down at r ~ \/6 as expected 

B.2.2 Intermediate field region r ^ 1 and r — (li^ 
The equation (^£)"^«'^)20mia! ^ q gj^Q^j^j^ ^^y^q background ( |B.9[ ) is 

2 



d 

dy 



y 



d 

dy 







(B.35) 



The solutions must match the far-field solution of the previous subsection, and the near-field 
solution of the next subsection, but are subject to no intrinsic boundary regularity requirements. 
The solution to this equation is 



ci + 



C2 



where Ci and C2 are constants. It is easy to check that the matching of (?/) with 



(B.36) 



sets Ci = y/k. It follows on general grounds that matching with the (as yet undetermined) near 
field solution forces C2 to vanish. This is because, were C2 to be non-zero, it would match onto a 
near field solution of order (|) (see the next subsection for details), violating the requirement 
that the solution has a smooth ^ — ?■ limit. Thus 



0r/^(y) = Vk 

This completes our intermediate field calculations at 0{\/9). 



(B.37) 
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At higher orders in 6, i.e. at 0{\/6 ^ ), the equation (B.35) generahses to 

2 



d 

dy 



d 



imid 



(y) 



-ymid 
(2m+l)/2 



(B.38) 



where P(tf+i)/2 



{y) is a source term determined by the perturbative procedure 
It turns out that the expressions ^'^^'^ 
?/ 1. At large y, (p'^2m+i)/2(y) scales like y 
an expansion in Oy"^. Further at ?/ ~ 1 



(2m+i)/2^y) increasingly singular as ?/ — )■ oo and 
2m other words the expansion of 0'"*'^ is really 



Lmid 



the 



,iy2hj) scales like ^y_i)m-i (m > 0), i.e. 
expansion for 0'"*°' is really an expansion in ^j^. The perturbation theory thus breaks down 
when y — 1 9 and y ^ \- 



B.2.3 Near field region r — or ?/ — 1^1 



We plug in (B.25) to the linearised equation 
The equation reduces to 



d 
dz 



d 



calculated in background (B.17). 



(B.39) 



/2 



Z) IS 



The near field solutions are subject to the requirement of regularity at the horizon z = 0. The 
leading order solution (after imposing regularity at the horizon) to this equation for (p" 

a constant. Matching determines the value of the constant to be y/k. 

<PTJz) = Vk 



(B.40) 



At OiVe ^^ ), the near-field scalar equation takes the form 



d 

dz 



d 



2Cz{l + 2z)-<P^^^^,y2{z) 



(2m+l)/2 



(B.41) 



As usual P(2m+i)/2^^) ^s ^ source term whose form is determined from the results of perturbation 
theory at lower orders. 



B.2.4 Summary 

We started this section with a small RN AdS black hole characterized by two parameters, a and 
a. In this subsection we have constructed an 0{-\f6) static correction to the starting vacuum 
solution. The correction lies entirely in the scalar field, and its amplitude is characterized by a 
new parameter k. 

The reader may be puzzled that we appear to have constructed a 3 rather than two param- 
eter set (as expected on general grounds) of black hole solutions. The resolution to this puzzle 
is that a will turn out not to be a free parameter, but will be determined as a function of a 
and k (this determination will, infact, justify the non interacting thermodynamic model of the 
previous section). However the constraint that determines a as a function of a and k shows up 
only at order 6'^/^, as we will see below. 
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B.3 Perturbation theory at 0{9) 

We now briefly outline the procedure used to evaluate the solution at 0{9). We proceed in 
close imitation to the previous subsection. The main difference is that at this order (and all 
integer orders) in the 6 expansion, perturbation theory serves to determine the corrections to 
the functions /, g, and A rather than (j). The procedure described here applies with minor 
modifications, to the perturbative construction at 0{6"^) for all m. The non-trivial equations 



at this order are the Einstein Equations (5.30) and the Maxwell equation (5.29) 



B.3.1 Far field region r ^ VO 



Plugging in (B.6) in (5.29) and (5.30), we get a set of three equations which we solve subject 



to boundary conditions (5.31). The equations take are 

32kr^ 



d 
dr 



r2(l + r2)2^°"*(r) 



dr 



3(l + r2)4 
2(1 + 2r2) 



32kr^ 



3(l + r2)5 



d 
dr 



^dA^\r 
dr 



32Jikr 



[l+r 



2^5 



(B.42) 



One of the integration constants in the first equation is fixed by the requirement that 



out/ 



is normalizable (see (5.31)). Another constant from the last equation is fixed by requiring that 
the gauge field y4^*(r) at infinity (see footnote 21). The remaining two integration constants 



(one in the first equation and one in the last) will be fixed by matching with the intermediate 
field solution below. With these conditions the solutions at 0{9) are 



^f*(r) 



out I 



'4\/6-9ci (r^ + 1 
18r2 (r2 + 1 




(B.43) 



where Ci and C2 are constants to be obtained by matching with the intermediate field solutions. 
To facilitate this matching, we note here the small r expansions of ^4°"*, /°"* and g""^^ 



A 



out 



out 



{9^ 



f 

out\ — l 



'l + 0{r)) + 9 
[l + Oir)) + e 



18r2 
-2a — C2k + 



8fc 
9 



8k 



-2a - C2fc + f 



2^6 

+ 0{r)]+0 (9^) 



(B.44) 



+ o{r) +o{e^) 



This completes the far-field calculations at 0{6). 
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At each higher integer order in 9, the far-field equations take the form (generahsations of 



(B.42) to 



d_ 
dr 



dr \ 1 + 
^f,dA^ 
dr V dr 



2(1 + 2r2 



-C\r) + Pi{r) 



(B.45) 



where P^, and are source terms determined by the solutions at all previous orders in 
perturbation theory. 

B.3.2 Intermediate field region r ^ 1 and r — (R'^ 



We plug in (B.14) to (5.29) and (5.30) to obtain the intermediate field equations. It turns out 



that the equations are slightly simpler when rewritten in terms of a new function 



Kmiy) = Voiy)9Tiy) + 



Vo{y) 



where 

Vo{y) = ( 1 

In terms of this function the final set of equations are 

dK,{y) 



dy 







2 d f 3c//l™^(?/)\ _ dKi{y) 



3 dy 



dy 



dy 



(B.46) 



dy 



{y^fr\y))='^yKM~2 



dA'^''^{y) 
dy 



These equations are all easily solved by integration, upto four undetermined integration con- 
stants (one each from the first and third equation, and two for the second). It will turn out 
that two of these constants are determined by matching with the far field solution while the 
other two are determined by matching with the near field solution. The solution at leading 
order is 



mid 



k bo- 



2y' 



fl 



mid 



9i 



mid 



+ 1 

r y 

%4 (V66i + 32/2(263 + &4)- 363) 



(B.47) 



3(2/2 
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Matching can be done by substituting y 
are 



and expanding around large r. These expansions 



A 



mid 



+ e 



2 V 2 



+ O (l/r^) 



jmid^ [r2 + l + 0(l/r3)] +e 
Matching determines 

Cl 



2a 



hk-- + (l/r^) 



(B.48) 



2a 



+ O (l/r^) 



4 2 



3 V 3 



(B.49) 



9 



To facihtate matching with the near field region in the next subsection, we substitute y = l+aOz 
and expand A™'^, Z™'^ and around z = 0. 

k 



^mid _ ^ I I 4 



(9 



+ 9^ (^-^a (Saoa - 6a + 2V6bik + 364A;) z + Aa^z^ + O (z^)^ + O (9^) 
0{9^) 



mid\ — l 



This completes our intermediate field calculations at 0{9) 



At all higher orders in 9, i.e. at 0{9"^^), the equations (B.46) generalise to 

dKmiv) 



dy 

2 d f.dAzHy) 



3 dy 



dy 



Si{y) + 



dKrn{y) 
dy 



(B.51) 



i^{y^f-\y))-sLiy)^2yK.M 



where S^{y), S^{y) and Sl^{y) are source terms determined by the perturbative procedure 
B.3.3 Near field region r — R<^R or y — \<^l 



To obtain the equations in the near field region, we plugin (B.25 ) in (5.29 ) and (5.30 ). However, 



as explained in ^:B.1.3, two linear combinations of these three equations (B.29) and (B.30) are 
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trivially zero. At leading order 0{6), we can therefore obtain only the relation for g^{^{z) in 
terms of A^^{z) and fl"'{z). This is 



2 1 d 



friz) 



(B.52) 



The functions /{"(-z) and A'-^{z) are obtained at which in turn determines g^{^{z). For 

the sake of completeness however, we list here the solutions for /["'{z), A'^i^z) and (7i"(-2) that 
are determined at 0{9^). 



friz) 



gTiz) 



4^/|A;((-42a + 8k + 21)z - 42alog(2;z + 1)) 
21a 



16A;(-42a + 8A; + 21)(2;z(856A; + 315(-42a + 8A; + 21)z) - 6615a{Az + 1) log(2z + 1)) 

416745a2 

28a2A;((-39690a + 808fc + 6615)^ + 6615a{4z + 1) log(2z + 1)) 



15(42a- 8^-21)3^2(2^ + 1)2 



(B.53) 



It can be verified that the (B.53) indeed satisfies (B.52). Matching can be done expanding by 



(|B.53|) at large z and comparing with (|B.50|). Expanding around large z, we get 

3 



6 ( ^Qaz + O [l/z) ]+e''\ -3x1 ^a'^z^ + ^^1^(7 - 16k)z 



+ 



= {Aah^ + 2o?{2 - ao)^ + O (1/z)) + {-\2o?z'^ + a^{-2a ~ 8k + 7)z'^+ 
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/clog 



42a - 8A; - 21 



3 V 42a^ 
16(210 - 821a)A;a 
^ 2205 
+0{l/z)) + 0{e') 



1930736A;2a 
1528065 



16 



z + ^a(42a -8k- 21)k 



2(3(a- l)a + l)a 

42a - 8A; - 21 



(log (- 



42az 



- 2 



(^*")-i =02 ^4^2^2 ^ Aa(42a - 8A; - 21)^ + O (l/z)^ + O {9^) 



(B.54) 



Matching clearly determines 



bi 

&4 



16 

y 



(B.55) 



The generalised near-field equations at every higher integer order in 9 were listed in (B.31) 



In the previous subsections, we have completely determined the first order corrections to 
the starting RN AdS black hole. However, in order to completely characterize the hairy black 
hole solution, we must also determine the leading order term in a as a function of a and k. 
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B.3.4 Determination of ao from the scalar equation at 0{9^^'^) 

In this section, we briefly present the far-fleld calculation at 0{9^^'^) which determines the 



leading order term in a as a function of a and k. We plug (B.6), (B.7) and (B.33) into the 
scalar equation 1)^0""* = and obtain an equation of the form (B.34). The source term in 



(B.34) is determined completely by the far-fleld solutions at 0{\/d) and 0{6), i.e. by known 
functions (p^^Kr), /°"*(r), gl^^{r) and y4°"*(r) and is explicitly 

IQy/kr 



tout I 



9 (r2 + 1 



-g [r^ (9r^ + 4(4A; + 9)r*^ + (44A; + 54)r^ + I2{k + 3)r2 + 9) 



(B.56) 



-9a (r^ + 1)^ (ttor^ + (ao - A)r^ + l) 



The solution to (B.34) is then given by 
^/k 



^3/2 ( 



189(1 + r2)2 
2k (-36rS 



9 log 



1 (-21a(ao - 4) + 8fc + 21) + 63 



1 — aao 



12a 
r2 + 1 



+ 



42r6 + eSr^ + 65r2 + 12) 



r2 (r2 + 1)' 



The small r expansion of 



Lout ( 



out I 



at C(03/2) 

^fke^/^{-21aaQ + Sk + 2l) 



+ O (rO) 



(B.57) 



(B.58) 



^3/21' J- g3^2 
Note that term of form (B.58) matches to a term of the form ^ in the intermediate fleld region. 



However, we know that the intermediate fleld solution at 0{\f9) is simply (f)'^^^{y) = \/k and in 
particular, possesses no term of the form ^ at large y. It follows that the corresponding term 



in (B.58) must vanish, giving 



1 8k 
a 21a 



(B.59) 



This determines the leading order term in a as a function of a and k and hence, completely 
determines our hairy black hole solution to 0{6) 

B.3.5 Summary 

To leading order the metric and gauge fleld of the hairy black hole is simply that of an RN AdS 
black hole parameterized by a and a. In this subsection we have completely constructed the 
leading order {0{6)) correction to the metric and gauge fleld of the hairy black hole solution. 
Importantly, we have also used one equation at 0{9^^'^) to determine a as a function of a and 
k (recall k parametrizes the scalar expectation value). In particular the results of this section 
establish that hairy black holes appear in a 2 parameter class, and justify the thermodynamical 
model of the previous section. We have carried out the perturbative construction of our solution 
to higher orders. In Appendix B.4 we note the full solution up to 0(6'^), and the far fleld solution 
for the metric and gauge fleld at 0{9^). 
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B.4 Listing of the final perturbative results 
B.4.1 Far field solution 

Scalar Field 



(r2 + 1)^ 



\/k (756a (r^ + 1)^ log + l) - 2 (378a (r^ + lf + k {A8r^ + 90r^ + Ar'^ - 17))) 



189 (r2 + 1)' 



8Vka{21a + Ak) ^, 

21 (r2 + if V 21(r2 + lf 



128aA;3/2 / 1 

log (r) + -Li2(-r' 



+ 



2aVk 



-378a (5r2 + 9) (r^ + l)^ 



189(1 + r2)6 

+4A; (2 (r^ + 2) (6r^ + 48r^ + 112r^ + 56r^ + 7l) + 107) + 567 (r^ + l) 
2y/k 



68762925r2 (r2 + 1) 



Jo -137525850a2 (4r^ + 8r2 + l) (r^ + l) 

+3465ar2 (r^ + l)'^ {2k (2520r^2 + 26460r^° + 104124r^ + 156336r^ + 91024r^ 
+29891r2 + 252071^ (r^ + + 1789^ + 19845 (r^ + l)^) + kr^ {k (2676720r^^ 
+10296810r^^ + 21833780r^° + 47080525r^ + 54295556r^ + 10445188r^ 
-8980072r^ - 5013164) - 5821200 (r^ + l)^ (9r^ + 21r^ + 13) 

Metric and Gauge Field 

+ 3r2 + 3) 



(B.60) 



Ar*(r) = 



2k Jl^^^ 



3(r2 + l)'^ 



3 A; 



105840ar2 (r^ + 3r^ + 3) (r 



l)'logf 4 + 1 



19845r2 (r2 + 1)^ 

+4410a (I7r^ + 44r^ + 18r^ - 36r^ - 9) (r^ + l)^ + (8A; (l31r^^ + 917r^° + 2751r^ 
+4270r^ + 4333r^ + 2541r^ - 747) - 6615 (r^ + l)^ + 3r^ + 3)) 

(B.61) 
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SkJla 



2205a (27r^° + 113r^ + 230r'^ + 306r^ + 225r2 + 27) (r^ + l)^ 



19845r2 (r2 + 1)^ 

(8A; (491 Ir^^ + 30825r^^ + 79527r^° + 106533r^ + 75845r^ + 22715/ 
-1995r^ + 735) + 6615 (3r^ + 7/ + 3r^ - 3) (r^ + 1)^)1 log (1 + ^ 



64 J lak {21ar^ (r^ + Sr'^ + 3) - 2A;) 



63r2 (r2 + 1)^ 



log' 1 + 



1_ 



-Lio 



63r2 (r2 + 1)^ V2 
k 

iT 



-r ) + log r + 



508295541600a^ (r^ + l)^Mog (ae'^ (2a 



8A; 
21 



31768471350y6r2 {r^ + 1) 
-2353220100a2 (221ri° + 901r^ + 1922r^ + 3105r^ + 2835r2 + 540) (r^ + l)^ 

+91476a (r^ + l)^ (sk (r'^ (^-22355646r2 + 588007r2 (r^ + 3r^ + 3) (r^ + l)^ - 
/ (2238406r^° + 16225148r^ + 50527008r^ + 87899161r^ + 93015307r^ + 60271386) 

-4384217) - 26145) - 77175 (r^ + l)^ (l42r^ + 342r^ + 222r2 + 9)) 

{-WP (((((3 (((55861681 (r^ + llr^ + 55) + 9207422235) 
+18541209786) + 26393430910) + 81319155082) + 61776138545) 
+31707662899) + 6259111287) + 748943493) 

-268939440A; ((((361 (r^ + 7r^ + 2l) + 12740) + 12677) + 7049) 
+1757) (r^ + 1) V 5294745225r2 (r^^ + llr^^ + 55r^^ + 164r^2 + 322r^° + 434r^ 
+406r^ + 260r^ + lOgr^ + 27)) + 15884235675r^] 

(B.62) 



/r(r) 



8k (r^ + 3r2 + 3) 
9 (r2 + 1)-^ 

64aA; (r^ + 3r^ + 3) log ^ (17^8 ^ 41^6 ^ 5^4 _ 54^2 _ ^g) 

9 (r2 + 1)=^ 27r2 (r2 + 1)^ 

_ 16A;2 (577ri2 ^ 3304^10 + 6972r« + 3248r6 - 7966r^ - 12966r2 - 8787) 

59535 (r2 + if 



(B.63) 
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32ak 

f?^Hr) = . [22050 (36r^2 + 180r^° + 365r^ + 434r^ + 402r^ + 270r^+ 

^ ^ ^ 59535r2 (r2 + 1)^ ^ ^ 

45) (r^ + iy + r^ (Sfc (((r^ (3 (((501 (r^ + 8) + 12949) + 21168) 
+17171) + 14420) - 2800) + 2730) + 6405) 

+19845 (r^ + 2) (r^ + r - l) (r^ + r + l) (r^ + l)^) log (^1 + 

256aA; (21ar2 (r^ + 3r2 + 3) - 2k) 



+ 



log^ 1 + ^ 

2048aA;2 ^1 ^ 2^ , 

Li2(-r ) + log r 

A: 



189r2 (r2 + 1) 

8k 



47652707025r4 (r^ + 1)^^ 



63536942700aV^ (r^ + l) log ( a^=^ ( 2a 



-5883050250^ (432r^^ + 2176r^^ + 4337r^° + 4531r^ + 3162r^ + 1818r^ 

+405r2 + 27) (r^ + l)^ - 22869ar2 (r^ + l)^ {4k ((l5536824r^ 
-588007r^ (r^ + 3r^ + 3) (r^ + l)^ + (420840r^2 + 4514806r^° + 19830323r^ 
+47657218r^ + 69544482r^ + 64721769r^ + 39221561) + 4359566) + 262605) 

+231525r2 (l2r^ + 67r^ + 134r^ + 113r^ + 3l) (r^ + l)^) 

+2kr^ {-k{{{{{3 (67989951r^° + 716821501r^ + 3397699745r6 + 9447771883/ 
+ 17302667096r2 + 22075104832) + 59792970787) + 39889757930) 
+ 19763824564) + 5824785222) + 2930284173) - 
800415 (6767r^^ + 48104r^° + 147252r^ + 256816r^ 

+271348r^ + 163008r^ + 42426) (r^ + l)^) 

(B.64) 



5 



' ^ ^ 9 (r2 + 1) 



52920ar2 (r^ + 3) (r^ + l)^ log + 1^ + 2205a (l7r'^ + 41r^ 



~ 59535 (r2 + l)'^ 

+18r2 + 54) (r=^ + l)^ + 2kr^ (577r^° + 3304r^ + 9912r^ + 29960r^ + 41930r2 + 7350) 

(B.65) 
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gl-\r) = — ^ [-2205a (r^ + l)' (5r« + 74r« + 108r^ - 72r2 + 9) + 

^ ^ ' 59535r2 {r^ + lfV \ ) \ ) 

8kr^ (3237r^^ + 18144r^° + 43932r^ + 70840r^ + 68950r^ + 16590r^ + 1260) 

+19845 (r' + 3) (r' + r)^] log (^1 + 

256aA; (21ar^ (r^ + 3) - 2k {3r^ + 1)) 2 A 1 
^ 189r2 (r2 + 1)^ °^ V ^ ^ 



2048A;2 (3ar2 + a) /U. , 2n , 2 



H ^ -Li2(-r^) +log^r 



8A; 

+ 



-635369427000^ (r^ + l)^^ log (^aO^ ~ ^ ~ -^)) 



47652707025r2 (r^ + l)^^ 
-588305025a^ (200r^° + 991r^ + 1403r^ + 936r^ + 1269r^ - 135) (r^ + l)^ 
-22869a (r^ + l)^ {4k ((588007r^ - 937666) r^^ + (4704007r^ - 6433583) r^^+ 
14 (I050007r2 - 1435547) r^^ + 672 (35007r^ - 56243) + 105 (l96007r^ - 400069) 
+175 (53 767r2 - 150641) + 1960 (907r2 - 6533) / - 1558830r^ - 21210) 

-231525 (r^ + l)^ (l3r^ + 40r^ + 38r^ + 24r2 + 3)) + 

2kr^ {k{{{{3 (67989951r^° + 716821501r^ + 3671382385r^ + 12206966475rV 
27226931600r^ + 45360308224) + 192955373842) + 174685548422) 
+52823111110) + 6712351338) + 800415 (6767r^° + 48104r^ + 147252r^ 

+229600r^ + 180460r2 + 58800) (r^ + l)^) 



(B.66) 



B.4.2 Intermediate field solution 

Scalar Field 



C/f (1/) = 

02'{y) = (378a log (a {y' - l) 9) + 189a {y' + 2) - 17^) 

\/k 

*^/"<»> = 68762925 fa^-l) I'^^^^^^^^"' + '"^^ + " 

+727650a ((y^ - l) (4158a - 428A; - 567) log(ae) + log (y^ - l) (y^(4158a - 428A; 
-567) + 1512a {y^ - l) log (a (y^ ^ l) ^) - 4914a + 572A; + 945) 

+18(42a - 8A; - 21) log (^^ (^2a - ^ - 1^ ^ + 756a {y^ - l) log2(a^) 

+756a {y' - l)'log (a (y^ - l) ^^) - 756a {y^ - l) log^ (^/ - l)) 

-6930a {2k (5565y^ - 3776y2 + 25207r^ (y^ - l) - 16909) + 19845 {y^ - 5)) 
+104A;(96407/c + 1455300) {y^ - l)] 

(B.67) 
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Metric and Gauge Field 



y2 



Ik 



8820a {2 log (^2a - ^ - 1^ + V log (a (7/^ - l) ^) - 41og(a) + 



fr\y) 



-6 log (y" - 1) - 21og(^) - 1) + (664A; + 2205) (y' - l)] 

16A; - 1) 



-13230a log (2« - ^ - 1) - log (a (y^ - l) ^) + 2 log( 



198451/4 

3 log {y^ - 1) + log(^)) + 6615a (y^ - l) + A; (292V - 2073)] 
8A;y4 (6615(a {y^ - 2) (y^ + 5) + 1) + 808A; - 26460a log {{y^ - 1) ^ (2a - |f - l))) 



B.4.3 Near field solution 

Scalar Field 



19845 (y2 - 1)'^ 



(B.68) 



C/2(^) 



^ (^189a (^2 log Q^a{2z + l)e\A2a -8k- 21)^ + log2(2;2 + 1)^ (B.69) 



3/2v-/ -^gg 

+378aLi2(-2^) + 63 (9 + tt^) a - 17A;) 
Metric and Gauge Field 



4W|A;(z(-42a + 8A; + 21) - 42alog(2^ + 1)) 
■^"('^ = Wa 

■fJ"'") = - '^^'tiwts!^ P-(315-(-12<' + 8* + 21) + 8564) (b.70) 

-6615(4a2 + a) \og{2z + 1)] 

in( _ 28a'^k{z{-39690a + 808fc + 6615) + 6615(4a^ + a) log(2^ + 1)) 
~ 15z2(2z + l)2(42a-8A;-21)3 

Matching also sets the parameter a as 



a 



\ 8k \ / 89224A;2 16A; 1 6568A; \, 

- + + h3a + -H 3]e + 0(e^) (B.71) 

Va 21a y V 1528065a 21a a 2205 / 



Near-horizon expansions of solution In order to see the breakdown of perturbation theory 
at extremality, we list the near horizon expansions of the near field solutions (in terms of the 
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AdS radial coordinate r) 
(f)^/^{r) ^-378alog ^^a^^(42a -8k- 21)^ - 63 (9 + tt^) a + 17k 

1764 v^A? (r- \^^) ^ 82320A; (r- ^ 



+ ^ ^ + ol(r- V^) 



^3(42a-8A; + 21)2 ^4(43^ _ gA; + 21)3^^ 



A.fik(-ma + Sk + 21)(r- VaO) 2352y6A; fr - VaO) 



2 



(a^)3/2(42a - 8A; + 21) 0^2(42^ _ gA; + 21)2 

65856 ( V6A;) (r - V^) ^ / . . 4^ 



Q3/2^9/2(42a _ 8A; + 21)3 



32A;(856A; - 6615a) (r - VaO) 32A;(-105a + 8A; + 21) (r - y/aO) ' 



19845aV2^3/2 3a3(42a - 8A; - 21)^3 

12544fc (r - VaO^ 

^ a5/2(-42a + 8k + 21)2^9/2 ^ ^ 



^m^^^ _4afc(-26460a + 808A; + 6615)(a^)='/2 ^ 56a2(6615(lla - 2) - 1616A;)A; 
45(-42a + 8A; + 21)2 (^r - 15(42a - 8fc - 21)3 

2352 ((735(58a - 9) - 808A;)A;(a^)=^/2) (r - VaO) / 
5((-42« + 8^ + 21W '- + 0{^{r-V^) 
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C Perturbative expansion for small charge black holes 



when = 



32 



C.l Far field solution 

Scalar Field 



(l + r2)2 



03/2(r) = 



2\/A? (r^ + 1)^ + 378a (r^ + 1)^ log (^^) + k (48r6 + 90r^ + - 17)) 



189 (r^ + 1)^ 
- 137525850a^ (4r^ + 8r2 + l) (r^ + l)* 



68762925r2(l + r2)io 
+ 104781600aA;r2 (r^ + l)^ Lis (-r^) + 6930aA;r2 (^25 207r^ (r^ + l)^ 

(4 (630r^ + 6615r^ + 26031r^ + 39084r^ + 22756) + 29891) + 1789) (r^ + l)^ 
- 727650ar^ (r^ + l)^ (^og (^^ + 1^ (^18 (r^ + l) ^ (21a + ik) log (^^ + 1^ 

-189a (5r^ + 9) (r^ + l)^ + 4A;r^ (r^ + 2) (6r^ + 48r^ + 112r^ + 56r^ + 7l) + 214A;) 
-288A; (r^ + l)^log2(r)) + /cV^ (r^ (r^ ((5 (535344r'' + 2059362r^ + 4366756r2 

+9416105) + 54295556) + 10445188) - 8980072) - 5013164) 



(C.l) 



Metric and Gauge Field 

2.fik (r^ + 3r2 + 3) 
Ar*(r) = ^ 

A-*(r) = 



3(r2 + l)' 
2k [2 
19845r2(l + r2)7V 3 



52920ar^ (r^ + 3r^ + 3) (r^ + l)^ log 



r2 + l 



(C.2) 



- 2205a (I7r^ + 44r^ + 18r^ - 36r2 - 9) (r^ + l)^ 

- Akr^ (r^ (((131 (r^ + 7r^ + 2l) + 4270) + 4333) + 2541) - 747) 



8A; (r^ + 3r2 + 3) 



9(r2 + l)" 
8A; 



52920a (r^ + 3r2 + 3) (r^ + l)^ log 



59525(1 + r2)7 
2205a (17r^ + 41r6 + 6r^ - 54r2 - 18) (r^ + 1)^ 



r2 + l 



(C.3) 



2k (577r^2 + 3304r^° + 6972r^ + 3248r^ - 7966r^ - 12966r2 - 8787) 
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grir) 



8kr^ (r2 + 3) 
8k 



- 52920ar^ (r^ + 3) (r^ + l)^log — 
\ / \ / y^,^ _|_ I 



59525(1 + r2)9 
+ 2205a (I7r^ + 41r^ + 18r^ + 54) (r^ + l)'' 

+ 2kr^ (577r^° + 3304r^ + 9912r^ + 29960r^ + 41930r2 + 7350) 

C.2 Intermediate field solution 

Scalar Field 

C/2 (y) = ^^ 

= (-378a log (a (y^ - l) ^) - 189a {y' + 2) + 17A;) 



(C.4) 



(C.5) 



Metric and Gauge Field 



2Jlk{y'-l) 



2205y2 



2205a (2 log r2a - ^) + V log (a (y^ - l) ^) - 2 log(a^e) + ^^'^^ 



/r''(z/) 



-6 log - 1) - 1) + 166A; (y' - l) 
Sk {y' - if 



3y4 



= 19^ [^22^°^ ((^' - ^) - I) + (2.^ log 



(a^) 



19845|/4 

+ (2y2 - 5) log {y^ - 1) - 3 log e) + (2 - 4^^) log a + 3 log {y^ - l) + log(^)) 
+ (y^ - 1) (6615a {y^ - l) + k (2929y^ - 2073)) 



Sky^ (^-26460alog (fg^zfy) + 6615a (y^ - 2) (y^ + 5) + 808A;) 



19845 (y2 - 1)"^ 



(C.7) 



C.3 Near field solution 

Scalar Field 



C/2(^) 



^ (378aLi2(-2^) + 189a(41og(a^) + log(2^ + l)(log(2^ + 1) + 2)) (C.8) 



189 

+63 (9 + TT^) a - 17k) 
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Metric and Gauge Field 

8^k{21az + 21alog(2z + 1) - 4:kz) 



21a 



friz) = [6615a(4(z - l)z - 1) + 6615{4az + a) log(2^ + 1 



416745a2 
+ 4^(4(208 - 315^)2 + 315)] 

7a^k{-6615a{10z + 1) + 6615(4a2 + a) \og{2z + 1) + 4fc(1462z + 315)) 



(C.9) 



30z2(2z + l)2(21a-4A;)3 
Matching also sets the parameter a as 



8k ^ ( 89224A;2 6568fc 



21a V 1528065a 



(C.IO) 



Leading order thermodynamics in perturbation theory Once we have the solutions 
from the previous subsections, the evaluation of their thermodynamic charges and potentials is 



straightforward. Using formulas (5.35), we find 

1024aA; 9232^^ 



M 

Q 

T 
S 



37r 

TT /3 
4V2 



^ 2a + 



8k 
"9" 



+ 



8k\ ,2 /1024afc 4192^^ 



189 59535 
+ 



3 2 /2,^ 

+ -\ -kO + 
2 TV 3 



V 189 ■ 59535 
^2 (4584195a2 + 4551624afc + 33784A;2 



1018710^6 

y^(21a - Ak) ^ 03/2 (_4584195a2 - 12701304aA; + 965368P 

2lTiy/a 
Ti^aOf/^ 



+ 



30561307rv^ 



+ o {e^) 

+ o [e^) 
+ o {e'/^) 



(C.ll) 



We have verified that these quantities obey the first law of thermodynamics 

dM = TdS + ndQ 

It is quite easy to show, that in terms of the shifted and rescaled mass Am and rescaled charge 
q, these formulae reduce to 



9a + 4k 
1 = -^r^ + C(0) 



3^6 



(C.12) 



Am 



567 



(63a + 8A;) + C(0) 



Positivity and reality of k constrains Am in the ranges 

7q2 

OiQ) > Am > — + 0{e) 
15 



(C.13) 



This range is shown in the Fig. 29 
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Am 

Ir 




Figure 29: Range of allowed values of Am for hairy black holes with = y at leading 
order in small 9 (blue shaded). Hairy black holes exist for all charges g > 0. The upper limit, 
Am = denotes the onset of superradiant instabilities (red line). This is also the extremality 
line for pure RN AdS black holes. The lower limit Am = —j^q"^ corresponds to the extremal 

hairy black hole with a = 2 (blue line). The green line denotes the soliton, Am = This 
phase diagram is quite similar to Fig. [16} the only difference being that the hairy black hole 
now exist at all charges g > 0. 
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D Perturbative expansion for = ^{1 -\- 6) 



D.l Far field solution 

Scalar Field 



03/2(r) 



2\/fc (378a (r^ + 1)^ - 378a (r^ + 1)^ log + l) + A; {ASr^ + 90r^ + - 17)) 



189 (r2 + 1) 



6 



- -18^1^ (5^^ + 9) + + 567 (r^ + 1)^ 

-4A; (2 (r^ + 2) (6r^ + 48r^ + 112/ + 56r^ + 7l) + 107)) log (^1 + 

8ayfc(21a + 4A;) , 2 A 1\ 64aA;3/2 
+ \ , ,^2 log 1 + - 01 n , ^2^2 (Li2(-r') + 2 V r) 



21(r2 + l)2 ^ V 21(l + r2)2 

- 137525850a2 {Ar^ + 8r^ + l) (r^ + l)' 



68762925r2(l + r2)io 
+ 3465ar2 (r^ + l)^ (2k (^25207r^ (r^ + l)^ + (4 (630r*^ + 6615r^ + 26031r^ 

+39084r^ + 22756) + 29891) + 1789) - 19845 (r^ + l)^) 

+ kr'^ {k (r^ (r^ ((5 (535344r^ + 2059362r^ + 4366756r^ + 9416105) + 54295556) 
+10445188) - 8980072) - 5013164) + 5821200 (9/ + 21r^ + 13) (r^ + l)"') 



(D.l) 



Metric and Gauge Field 

2Jjk (r^ + 3r2 + 3) 



3 (r2 + 1)^ 
^k 



^r(r) -- 

AT\r) = -X^ ^ I" - 105840ar' + 3r^ + 3) (r' + l)^ log + 

^ 19845r2(l +r2)7 [ ^ M J ^ ^^.2 J p ^) 

- 4410a (I7r^ + 44r^ + 18r^ - 36r2 - 9) (r^ + l)^ 

- {Sk (r^ (((131 (r^ + 7r^ + 2l) + 4270) + 4333) + 254l) - 747) 
+6615 (r^ + 3r^ + 3) (r^ + l)^) 
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/r(r) 



8k (r^ + 3r^ + 3) 



9(r2 + l)^ 
8A; 



52920ar2 (r^ + Sr^ + 3) (r^ + l)^ log 



r2 + l 



grir) = 



59535r2(l + r2)7 
2205a (I7r^ + 41r^ + 6r^ - 54r^ - 18) (r^ + l)^ 

2kr^ (577r^2 + 3304r^'' + 6972r*^ + 3248r^ - 7966r^ - 12966r2 - 8787) 

8kr^ (r^ + 3) 
9(r2 + l)^ 

8k 



52920ar^ (r^ + 3) (r^ + l)^log (J^ + 1^ 



59535r2(l + r2)9 
+ 2205a (I7r^ + 41r^ + 18r2 + 54) (r^ + l)'' 

+ 2kr^ (577r^° + 3304r^ + 9912r^ + 29960r^ + 41930r2 + 7350) 

D.2 Intermediate field solution 

Scalar Field 



^1/2 

imid ^ 

3/2 - 



y/k (-378a log (a {y^ - l) ^) - 189a {y^ + 2) + 17A;) 



Metric and Gauge Field 



2Jlk{y'-l) 



2205^/2 



(y2 - 1) (8820a (y2 + 1) + 664A; + 2205) 



8A; 



(D.3) 



(D.4) 



fr\y) 



17640a y- log (^2a - — + 1 ) - 2|/2 log (a (i/^ - l) ^) + log(a^^) + 3 log (i/^ - l) 
8A; - 1)2 



3|/4 

8fc(y2 - 1) 



13230a (y^ - l) + A; {b8b8y'^ - 4146) - 6615 



198457/4 

26460a r31og (y^ - l) - log ('2a - ^ + 1 ) - 2y2log (a (y^ - l) 9) +2\og{a''e) 



gf\y) = 



gT'iy) 



8ky^ (-26460a log 



e(42a-8fc+21) 
21(2/2-1) 



6615a - 2) (7/2 + 5) + 808A;) 



19845 (7/2 - l)'^ 



(D.5) 
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D.3 Near field solution 

Scalar Field 

<f>T/,{z) = Vk 

0372(^) = -J^v^(l89a(^21og(^^a^'(2z + l)(42a-8fc + 21)^ + \og\2z + l)^ (D.6) 
+378aLi2(-2z) + 63 (9 + tt^) a - 17k) 
Metric and Gauge Field 

A.[\k{z{A2a - 8k + 21) + 42a log(2;z + 1)) 
= 

, = JK12«-Sk + 2l) rj2gj^2,42<, _ 8fc + 21) + 315(42a -Sk + 21) 

416745a^ 

+13230(4az + a) log(2z + 1) + 2(6615 - 1712A;)z] 

14a2fc(-13230a(6z - 1) + 13230(4a;z + a) log(22 + 1) + 8A;(2022 - 315) + 6615) 

15^2(2^ _^ l)2(42a -8k + 21)3 

(D.7) 



Matching also sets the parameter a as 

8A;-21 /fHi^ + ^ + l „ 6568A; 



a 



^gl 1528065 21 ^3a+^^+3 (D.8) 



21a \ a 2205 
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